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ABSTRACT 


An  evaluation  of  the  NIST  FDS  model  was  conducted  with  particular  attention  for  its  use 
in  predicting  flame  spread  on  surfaces.  Over  the  course  of  this  investigation  the 
computational  model  changed  from  combustion  depicted  by  particles  to  a  mixture  fraction 
based  combustion  model.  The  study  pertains  to  version  2.0  released  on  December  4, 
2001. 

Three  aspects  were  considered  in  the  study.  First,  we  studied  the  evaluation  of  the  code  to 
predict  a  combusting  plume.  Second,  the  code  was  applied  to  a  fire  plume  adjacent  to  a 
vertical  wall,  and  then  flame  spread  on  the  wall.  Third,  a  complementary  investigation  of 
an  improved  algorithm  for  convective  heat  transfer  at  a  surface  was  developed.  The  first 
two  studies  resulted  in  M.S.  theses.  Damian  Rouson  of  CCNY  performed  the  third  study. 
The  thesis  by  Ma  on  the  axi-symmetric  plume  was  previously  transmitted  and  will  not  be 
included  here.  However,  a  recently  accepted  paper,  based  on  the  thesis  with  updated 
results  is  included. 

The  general  conclusions  are  that  the  FDS  code  is  very  good  for  computing  the  fluid 
dynamics,  entrainment  and  flame  height.  The  temperature  in  the  combustion  region 
appears  to  be  over-estimated  at  the  base  of  the  geometry  considered,  and  any  related  heat 
flux  is  consequently  over-predicted.  The  temperature  results  are  grid  dependent.  A 
computation  of  flame  spread  on  vertical  PMMA  gave  mixed  results. 

The  code  was  benchmarked  against  fire  plume  correlations  after  a  review  of  the  literature 
to  obtain  the  most  general  results.  Most  of  the  experimental  correlations  have  some 
deficiencies,  and  should  be  improved.  Particular  attention  needs  to  be  given  to 
temperature  measurements  in  the  flame  since  these  are  generally  under-estimated  due  to 
radiation  error.  The  wall  heat  flux  and  flame  spread  comparisons  were  made  against  data 
we  viewed  as  quality  data. 

The  algorithm  developed  by  Rouson  is  based  on  the  theoretical  formulation  by  Howard 
Baum,  and  has  not  been  tested  in  the  FDS  code. 
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Executive  Summary 

The  objective  of  this  study  was  to  investigate  the  ability  to  compute  flame  spread  on 
walls  using  a  field  model.  The  model  selected  was  the  Large  Eddy  Simulation  (LES) 
developed  by  NIST/BFRL  that  is  entitled  the  Fire  Dynamics  simulator  (FDS). 
Version  2.0  was  used  in  the  final  analysis  presented  in  this  report;  however,  this  work 
was  originally  begun  with  version  1.0.  In  version  2.0,  the  combustion  is  computed 
using  a  mixture  fraction  formulation,  while  in  version  1.0  the  combustion  was 
specified  by  particles  that  released  energy  until  their  burnout.  The  time  release 
depended  on  the  nature  of  the  fire  being  simulated.  Thus,  version  1.0  was  not  suited 
for  predicting  fires  of  an  arbitrary  nature  as  those  resulting  from  a  particular  geometry 
or  spread. 

Particle  Energy  Release  Rate.  Version  1 .0,  with  the  prescribed  particle  energy 
release,  is  still  a  useful  approach  when  information  about  the  fire  is  available.  The 
algorithm  for  describing  the  time-release  of  the  particle  energy  is  based  on  the 
McCaffrey  fire  plume  correlation  for  axi-symmetric  fires,  and  a  constant  energy 
release  rate  applied  to  each  particle.  For  fires  of  other  geometries,  expressions  must 
be  available  for  representing  the  characteristic  velocity  and  flame  length,  in  the  least. 
In  a  combustion  system,  the  energy  release  rate  is  not  uniform,  but  controlled  by  the 
rate  of  oxygen  entrainment.  A  more  general  approach  to  selecting  the  energy  release 
rate  of  the  particles,  and  to  illustrate  its  dependence  on  the  fire  geometry  is  presented 
below  since  it  was  not  included  in  past  progress  reports. 

Represent  the  energy  release  rate  of  the  fire  in  terms  of  the  selected  particle  release 
rate  and  its  distributed  energy  release  rate  over  the  time  for  combustion. 

£?  =  ",.  /n  <iP  (t)dv  =hp  qp  (rj)dr] 
where  qp  is  the  particle  energy  release  rate, 
th  is  the  burning  time, 
zf  is  the  flame  length, 
hp  is  the  selected  particle  rate. 


The  velocity  of  the  particles  generally  depends  on  their  launch  site.  But  if  they  are 
based  on  a  characteristic  plume  velocity,  usually  taken  as  the  centerline  velocity,  then 
the  trajectory  can  be  based  on 


dz 

w  =  —  or  t 

dt 


=  f  —  where  w  would  be  known  as  a  function  of  z. 
w 


The  rate  of  entrainment  into  a  fire  plume  can  be  approximately  represented  in  terms 
of  an  entrainment  velocity  proportional  by  a  constant  (ajto  the  characteristic  velocity 


and  the  perimeter  of  the  fire  plume,  P(z).  It  follows  that  the  rate  of  mass  entrainment 
is 

me  =  £ap^wPdz. 

The  rate  of  combustion  energy  in  the  plume  depends  on  the  oxygen  that  burns,  and 
this  has  been  empirically  found  to  require  n  times  the  needed  stoichiometric 
requirement  with  n  generally  found  as  about  10.  The  combustion  rate  is  related  to  the 
entrainment  as 


0  =  7,  In. 


The  particle  energy  release  rate  must  follow  accordingly  as 


q  =  apJfoxxkhcw2P  l(nhp). 

"  dz  p 


The  velocity  in  the  combustion  region  of  a  fire  plume  is  found  to  be  proportional  to 
z"2.  The  distribution  of  P  with  z  depends  on  the  dimension  of  the  firebase  relative  to 
the  flame  height.  For  short  flames,  P  is  nearly  constant  and  related  to  the  base 
dimension;  and  for  tall  flames,  P  is  proportional  to  z-  Therefore,  the  particle  energy 
release  rate  should  vary  as  z  or  z 2-  For  flames  that  have  the  flow  constrained,  as  for 
example  corner  fires,  P  would  decrease  accordingly. 

While  this  result  for  the  particle  energy  rate  is  more  complex  than  in  version  1.0  of 
the  FDS ,  it  is  more  general  and  may  find  utility;  or  in  the  least,  shows  the  dependence 
of  the  particle  energy  rate  on  the  plume  geometry. 


Fire  plume  modeling.  The  objective  of  this  study  was  to  assess  the  application  of  a 
CFD  model  for  the  prediction  of  fire  spread.  The  model  selected  for  investigation 
was  the  NIST  FDS  code,  and  the  objective  became  to  establish  its  accuracy.  Before 
the  general  problem  of  fire  spread  could  be  examined,  it  was  felt  that  the  accuracy  of 
the  FDS  code  should  be  assessed  with  simpler,  but  related  problems.  Since  a  fire 
plume  is  usually  the  ignition  source  and  driving  heat  transfer  for  fire  spread  on 
adjacent  objects,  the  fire  plume  became  the  first  issue  for  examination.  Moreover, 
much  data  exist  for  axi-symmetric  fire  plumes.  These  data  became  our  focus  in  order 
to  assess  the  FDS  accuracy.  Our  approach  was  to  examine  the  array  of  data  and 
correlations  for  fire  plumes  and  establish  the  “best”  bases  for  comparison.  “Best” 
was  decided  on  the  most  comprehensive  fit  to  data,  and  on  the  completeness  of  the 
correlation  variables.  Often  both  criteria  could  not  be  satisfied  since  most 
correlations  do  not  include  all  of  the  relevant  variables.  The  details  are  found  in  the 
Section  1. 

For  example,  centerline  plume  temperature  above  the  flame,  in  general,  depends  on 


2 


where 


T-T 


=  f(Xr,zlz',Dlz*) 


D 


{p.C'T.'fej 

The  best  correlation  above  the  flame  was  found  based  on  the  work  of  Heskestad: 


T-Tx 

T 


=  cT{\-xry 


m  Z  -  Z„ 


-5/3 


where 


-^-  =  -1.024  +  1.37 


Z  7 

gives  the  virtual  origin,  and  Heskestad  gives  CT  =9.1 

In  the  flame  region,  the  centerline  temperature  is  approximately  constant  and  depends 
on  radiation  fraction  and  somewhat  on  the  fuel.  It  is  shown  that 


T-T, 


=  4.89(1  -  X,)  or  T  -  Tx  =  1455"C(1  -  Xr). 


This  suggests  that  the  maximum  temperature  in  turbulent  flames  with  no  radiation 
loss  could  be  as  much  as  1500"C. 

The  experimental  results  for  flame  height  are  shown  below  for  the  common 
correlation  of  Heskestad,  correlations  from  Zukoski  over  various  Q*  values  including 
values  below  0.1,  and  that  of  Quintiere  and  Grove  which  included  the  effect  of  flame 
radiation  fraction. 


Figure  1 .  Flame  height  correlations 
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The  FDS  computations  show  good  results  for  flame  height  as  depicted  below.  The 
computed  results  use  a  criterion  of  the  percent  of  combustion  completeness  to  define 
the  flame  height  in  the  computation.  The  computation  results  for  low  Q*  <  0.5 
deviate  markedly  from  the  correlation.  However,  the  data  in  this  region  are  scant,  and 
a  more  thorough  experimental  study  is  needed  to  properly  examine  the  computations 
here 


Figure  2.  Computational  results  for  flame  height. 


The  computed  results  for  the  centerline  Froude  number,  Fr  =  w2  h 


T-L 


\ 


gz,  are 


shown  in  Figure  3.  Inviscid  theory  gives  a  value  of  2,  and  experimental  values  of 
about  1.3- 1.9  are  generally  found  in  the  literature.  The  FDS  computation  does  a  good 
job  a  reproducing  these  values  and  shows  a  distinct  effect  of  Q*  that  should  be 
examined  experimentally 


Figure  3.  Fr  computed  on  the  centerline. 
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The  rate  of  entrainment  over  the  flame  height  was  also  computed,  and  the  entrainment 
rate,  normalized  with  the  stoichiometric  required  airflow  rate,  was  plotted  against  the 
normalized  height.  Figure  4  shows  that  for  <2*>0.5,  the  amount  of  excess  air 
entrained  roughly  10  +/-  2  as  is  generally  found  from  experiments.  This  is  a  measure 
of  the  “unmixedness”  of  the  fuel  and  the  air  in  the  fire  plume.  Lower  values  of  Q* 
exhibit  much  lower  excess  air  ratios.  Again,  insufficient  experimental  data  do  not 
allow  a  definitive  statement  on  the  FDS  accuracy  in  this  Q*  range.  However,  FDS 
appears  to  do  a  good  job  at  predicting  the  entrainment  rates,  suggesting  that  the  large 
eddy  structure  captured  in  the  model  may  be  most  significant  for  entrainment. 


Figure  4.  Normalized  entrainment  rate  in  the  flame 


Perhaps  the  weakest  ability  of  the  FDS  code  lies  in  the  ability  to  predict  the 
temperature  in  the  flame  region.  This  is  likely  due  to  the  inability  of  the  code  to 
resolve  the  sub-grid  combustion  phenomena.  The  computed  results  compared  to  the 
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The  computed  temperatures  tend  to  overestimate  the  correlation  values  in  the  flame 
zone  for  Q*<  1,  and  show  the  effect  of  introducing  “cool”  fuel  for  large  Q*  values. 
The  comparison  in  the  far  field  plume  region  is  very  good.  The  combustion 
temperatures  can  be  overestimated  by  about  200  degrees.  Again,  good  measurements 
in  this  region  are  lacking  since  radiation  effects  influence  the  measurement,  and  data 
are  lacking  for  low  Q*  fires  in  general. 


In  all  of  our  computations,  we  used 

R*  = 


max(&c,<5y,<5z) 


as  a  criterion  for  selecting  the  grid  size  to  insure  sufficient  cells  over  the  characteristic 
flame  length, 

Z  ={Q/pxcnT«J>)215. 

We  used  R*  =  0.05  or  20  cells  over  the  characteristic  flame  length.  In  FDS  2.0  the 
stoichiometric  mixture  fracture  is  adjusted  to  compensate  for  the  code  to  resolve 
combustion  is  large  grid  cells.  This  adjustment  is  given  as 


min  1  ,C- 


where  C  is  an  empirical  constant  to  be  used  for  all  fire  scenarios  (here  it  is  0.025  in 
FDS2.0),  is  the  grid  size  and  D*  (=  herein).  Thus,  for  all  R*  >  0.025  there  is  an 
adjustment;  and  our  criterion  just  fell  into  this  category.  An  experimental  set  of 
computations  for  a  fire  of  72.5  kW  with  varying  grid  cell  selections  was  done  to 
investigate  the  effect  of  grid  on  the  flame  temperature.  These  results  are  shown 
below.  The  temperature  in  the  flame  appears  to  increase,  roughly,  as  the  grid  size  is 
decreased.  Except  for  the  coarsest  grid,  the  far  field  plume  temperatures  appear  to  be 
grid  independent.  The  grid  dependence  is  not  understood,  and  suggests  a  weakness  in 
the  current  combustion  model. 


Table  3.  Simulation  result  of  grid-dependent  cases. 


No. 

Grid  number 

Domain  size 
(m) 

Grid  size 
(cm) 

Computed 
combustion 
rate  (kW) 

Total 

Radiative 

fraction 

Maximum 

Temp.(°C) 

Estimated 

Zf.cf/Zf 

1 

32x32x72 

1.5x1. 5x3.5 

4.69 

65.0 

0.38 

1205 

0.18 

2 

48x48x96 

1.5x1. 5x3. 5 

3.13 

64.5 

0.40 

1288 

0.28 

3 

64x64x144 

1.5x1. 5x3.5 

2.34 

66.6 

0.38 

1000 

Esam 

4 

32x32x72 

1 .5x1 .5x3.5 

Non-uniform 

74.9 

0.46 

1374 

5 

32x32x80 

,54x. 54x1.35 

1.69 

72.8 

0.41 

1219 

0.51 

6 

64x64x160 

.54x.54xl.35 

0.84 

86.4 

0.43 

1347 

1.00 

6 


Width  (m)  Height  (m) 


Figure  6.  Temperature  dependence  on  grid  size. 


Wall  fires.  The  next  task  was  to  investigate  the  accuracy  of  the  FDS  code  in 
predicting  the  heat  transfer  to  wall  subject  to  an  adjacent  pool  fire.  The  experiments 
of  Back  et  al.  served  as  an  experimental  basis.  The  experimental  arrangement  is 
shown  in  Figure  7.  The  predicted  heat  flux  was  in  fair  agreement  with  the  magnitude 


of  the  measured  values,  but  was 
distributed  too  tightly  as  shown  in 
Figure  8.  However,  the  flame 
height  was  reasonably  predicted 
according  to  Figure  9.  This  is 
attributed  to  the  over  Drediction 
of  the  temperature  in  the  region 
closest  to  the  burner.  The  same 
results  characterize  the  free 
plume  predictions.  It  appears 
where  more  fuel  exists  the  burn 
rate  is  too  high,  but  the 
distribution  of  burning  occurs 
over  the  proper  region.  This 
Square  propane  might  be  interpreted  that  burning 

burner  w/  Q=50~500  rate  m  the  cell  is  too  high  in  fuel 

kW  edge  0  28-0  7  m  rich  regions,  hut  the  coarse  scale 

mixing  of  the  fuel  and  air  is 
consistent  and  a  good  prediction 
of  flame  height  results. 

Propane  supply 


1.25  cm  thick 
gypsum  board  back 
by  steel  frame 


Thermocouple  1.25 

cm  r - surface; 

heat  flux  meter  on 
surface;  1 5  cm  apart. 


Figure  7.  Wall  heat  flux  measurements. 
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Flame  Spread.  The  flame  spread  computations  were  compared  to  vertical  spread  over 
a  5  m  slab  of  PMMA.  These  experiments  were  conducted  by  FMRC.  There  was 
some  difficulty  in  simulating  the  ignition  process,  and  the  modeling  of  this  was 
obverted.  Instead,  synchronizing  at  a  low  energy  rate  indicative  of  the  incipient  fire 
began  the  comparison.  The  experimental  arrangement  is  shown  in  Figure  10.  The 
comparisons  for  the  spread  rate  and  energy  release  rate  during  spread  are  shown  in 

Figures  11  and  12, 


To  Fire  Product  Collector 


1.2m 


Flame  Extension 
Board  (24  gauge 

steel) _ 

PMMA 


Marinite 


Flow  holder 
(24  gauge  steel 
0.6  m  wide) 


3  m 


5  m 


respectively. 
predictions  are  very  good. 
However,  a  comparison  at 
the  quasi-steady  burning  rate 
once  all  of  the  5  m  was 
burning  shows  higher 
burning  rates  than 
measurements  in  Figure  13. 
The  measurements  were  not 
taken  from  the  flame-spread 
experiments  since  none  were 
measured  there.  The 
measurements  were  taken 
from  a  1.2  m  tall  PMMA 
study  by  Orloff.  Estimates 
were  also  made  from  a 
correlation  developed  by 
Ahmad  and  Faeth.  When  a 
smaller  vertical  slab  of 
PMMA  was  modeled,  we 
believe  this  higher  heat  flux 
contributed  to  a  “runaway” 
causing  an  accelerating 
spread  rate,  and  termination 
of  the  run.  That  result  is 
shown  in  Figure  14  along 
with  a  prediction  from  the 
SOFIEcode. 


Load  Cell 


Figure  10.  Flame  spread  experiment 
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Figure  1 1 .  Flame  spread  rate  comparisons. 
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Figure  12.  Energy  release  rate  comparison. 
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Figure  13.  Burning  rate  comparison. 
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Figure  14.  Prediction  for  spread  on  1  m  vertical  PMMA 


Convection  heat  transfer  at  the  wall.  A  sub-task-  of  this  work  was  to  examine  a 
method  to  improve  the  convective  heat  transfer  at  the  wall.  Currently  the  code  uses  a 
standard  heat  transfer  correlation.  The  improved  algorithm  was  based  on  a  model  of 
Howard  Baum  that  solved  the  exact  unsteady  energy  equation  with  a  uniform  velocity 
parallel  to  the  wall  based  on  the  closest  grid  cell.  This  is  depicted  in  Figure  15. 
Damian  Rouson  examined  the  model,  and  an  algorithm  was  developed.  Comparisons 
were  also  made  with  exact  solutions. 
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Figure  15.  Convection  at  the  wall. 

Concluding  Remarks.  Details  of  the  findings  are  found  in  the  following  Sections. 
They  should  be  examined  for  a  complete  understanding  of  the  analyses.  In  general, 
our  experience  with  the  current  FDS  code  was  good  in  terms  of  its  ability  to  make 
predictions.  The  weak  point  appears  to  be  the  combustion  model  that  causes  higher 
temperatures  and  heat  fluxes  due  to  the  fire,  especially  in  the  fuel  rich  region.  The 
predictions  for  flame  height  were  good  suggesting  an  accurate  representation  of  the 
“unmixedness”  character  of  turbulent  fires.  The  thermal  and  fluid  dynamics  outside 
the  combustion  zone  appeared  to  yield  good  results  for  a  free  fire  plume.  While  one 
application  of  flame  spread  gave  good  results,  another  produced  an  incorrect 
“runaway”  condition.  However,  the  runaway  is  not  inconsistent  with  the  physics  of 
upward  flame  spread 
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Abstract: 

The  objective  of  this  work  has  been  to  determine  the  accuracy  and  limitations  of  a  new 
version  of  Fire  Dynamics  Simulator  (FDS),  developed  by  McGrattan  et  al.,  on  axi- 
symmetric  fire  plumes.  The  current  version  uses  FES  for  turbulence,  a  mixture-fraction- 
based  infinitely  fast  chemistry  model  for  combustion,  and  a  constant  radiative  loss 
fraction.  These  sub-models  have  been  tested  for  unconfined  fires  of  different  sizes,  based 
on  a  dimensionless  heat  release  rate  Q*d  in  the  range  of  0.1  to  10.0,  which  covers  most 
natural  fire  scenarios.  No  adjustment  of  constants  or  algorithms  in  the  model  FDS2.0 
have  been  made.  An  examination  of  plume  theory  is  made  first  to  find  the  benchmark 
correlations.  This  shows  a  generalization  for  a  collection  of  correlations  based  on  theory, 
and  which  might  be  “the  best”.  Using  the  characteristic  length  as  the  scaling  factor,  it  is 
found  that  the  optimum  resolution  of  a  pool  fire  simulation  is  around  0.05.  With  this 
resolution,  the  flame  height  prediction  is  found  to  fit  well  with  flame  height  correlations. 
Some  other  parameters  such  as  temperature  and  mixture  fraction  are  found  to  be  close  to 
the  empirical  estimations  at  flame  tips.  The  Froude  number,  which  describes  the  relative 
strength  of  momentum  and  buoyancy,  falls  within  the  measurement  range  of  many 
researchers.  The  simulation  also  reveals  that  the  temperature  near  the  burner  is  over¬ 
predicted,  while  the  centerline  temperature  and  velocity  in  the  non-combusting  region  is 
predicted  well. 

Keywords: 

Fire  plume  dynamics,  flame  height,  entrainment,  large  eddy  simulation 
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LIST  OF  ABBREVIATIONS 


Nomenclature 

c  specific  heat  of  air  (J/gK) 

C  comstant  in  FDs  to  adjust  the  stoichiometric  mixture  for  grid  size,  equation  4b. 

q  constant  for  temperature  distribution  in  plume  region 

C  constant  for  velocity  distribution  in  plume  region 

q  constant  for  maximum  temperature  in  flame  region 

D  effective  diameter  of  burner/  species  diffusivity 

D *  the  local  characteristic  length  with  reference  to  gridsize 

zf  flame  height  (m) 

I  radiation  intensity 

Ib  radiation  blackbody  intensity 

m"  fuel  mass  loss  rate  per  unit  area  (kg/s/m2) 

n  entrainment  number,  the  ratio  of  flame  entrained  air  to  stoichiometric  air  needed. 

N  dimensionless  group  for  fire  power,  equation  (3) 

Q  fire  power  (kW) 

n*  non-dimensional  heat  release  rate  based  on  burner  diameter 

1  heat  release  rate  per  unit  volume,  HRRPUV  (kW/nr) 

R*  resolution  of  simulation,  gridsize/z* 

T  temperature  ("C) 

t  time  (s) 

U  integrated  radiant  intensity  (kW/m2) 

u,v,w  directional  velocities  (m/s) 

x,y,z  spatial  coordinates 

s  stoichiometric  oxygen/fuel  mass  ratio 

sa  stoichiometric  air/fuel  mass  ratio 


Xr  Radiative  Loss  Fraction 

Yf  mass  fraction  of  fuel 

Ya  mass  fraction  of  oxygen 

z  vertical  coordinates 

z  t  flame  height 

Z  mixture  fraction 

Z  averaged  mixture  fraction  at  flame  tip 

7  stoichiometric  mixture  fraction  at  flame  sheet 

•y< 

7  effective  mixture  fraction  at  flame  tip 

z  the  characteristic  length  of  the  fire  with  reference  to  the  source  diameter 
Ah,  fuel  heat  of  combustion  (kJ/kg) 

Ah„  energy  released  per  unit  mass  oxygen  consumed 
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absorption  coefficient  (m'1) 
density  (kg/m3) 


1.  Introduction 

As  a  basic  phenomenon,  the  buoyant  plume  has  always  attracted  attention 
throughout  fire  research  history.  A  lot  of  investigators  have  contributed  to  the  fire 
plume  research  with  well-designed  experiments.  Among  them,  Morton  et  al,  m,  Yokoi 
|21,  Thomas  et  al.131,  McCaffrey  |41,  Zukoski  et  al.151,  Hasemi  et  al,161,  Cox  et  al.171, 
Heskestad  |S|,  Delichatsios  |91,  Kung  et  al.|l0!,  Hamins  et  al.1"1,  etc.,  made  significant 
contributions  to  the  plume  theory  and  its  understanding.  These  studies  represent  simple 
physical  Froude  modeling  and  dimensional  analysis  to  model  buoyant  turbulent  fires. 
Most  of  these  contributions  were  summarized  by  Beyler  il21,  McCaffrey  |l31,  Heskestad 
|Hi,  Delichatsios  [l5),  Zukoski |l61,  and  most  recently  Quintiere  etal.1'71. 

With  the  advancement  of  computer  technology  and  computational  methods, 
more  and  more  research  work  shifted  to  the  numerical  simulation  of  fire  plumes.  Most 
numerical  simulations  are  based  on  the  k  -  £  turbulence  model  (Reynolds  average 
simulations,  RAS).  Tamanini  1181  used  an  improved  version  of  the  k  -  e  -  g  model  of 
turbulence  to  simulate  buoyancy  controlled  turbulent  diffusion  flames.  His  combustion 
model  is  infinitely  fast  chemical  kinetics  and  local  radiant  emission  is  simplified  to  be  a 
constant  fraction  of  total  energy.  You  and  Faeth  1191  used  a  model  similar  to  Tamanini’s 
k  -  Er  g  model,  only  lessening  the  effect  of  buoyancy,  applied  to  flow  which 
buoyancy  and  low  Reynolds  number  effects  were  small.  They  provided  the  most 
complete  comparison  between  prediction  and  experiment  in  buoyant  flames.  Crauford 
et  al.  |20t  used  k  -  Emodel  and  flamelet  combustion  to  simulate  a  25  cm  turbulent 
natural  gas  flame.  Adiga  1211  used  the  k  -  Emodel  and  eddy-dissipation  concept  for 
combustion.  His  radiation  model  is  based  on  the  flux  model  and  temperature- weighted 
gray  gas  model.  Both  used  a  parabolic  flow  field  assumption. 

Beginning  in  the  late  1980’s,  some  general-purpose  commercial  CFD  codes 
were  available  in  the  market.  Some  of  them  were  modified  and  applied  in  fire  research, 
such  as  JASMINE,  KAMELEON,  FLUENT,  and  SOFIE,  etc.  It  is  pointed  out |22!,  the 
k  ~  e  models  tend  to  over-predict  velocities  and  temperatures  along  the  centerline  axis 
of  buoyant  plumes  and  consequently  under-predicts  the  plume  width.  Some  suggestions 
were  made,  such  using  as  the  Algebraic  Stress  Model  (ASM),  to  account  for  the 
influence  of  buoyancy.  Nam  and  Bill  1231  modified  the  standard  k  -  e  model  in 
PHOENICS  by  adjusting  its  constants  to  study  the  pure  thermal  plume  above  the  flame 
region.  No  combustion  model  was  used.  The  velocity  half-width  based  on  velocity  and 
temperature  was  then  compared  with  those  of  standard  k  -  e  model  and  experimental 
data  of  large  fires  by  Kung  et  al.  |in|.  This  study  1231  is  part  of  a  research  program  devoted 
to  numerical  simulation  of  the  interaction  of  fire  and  sprinkler  spray.  Sinai  et  al,  1241  used 
FLOW3D  to  model  the  unconfined  pool  fires  with  emphasis  on  the  flame  geometry 
above  a  round  20  m  kerosene  pool  fire  subjected  to  a  crosswind.  The  role  of  the  pool 
shape  and  ambient  turbulence  has  been  investigated.  His  model  is  based  on  the 
buoyancy-modified  k  -  e  turbulence  formulation,  the  Eddy  Break-up  combustion 
model  and  a  gray-medium  representation  of  thermal  radiation.  Miles  et  al,  1251  used 
JASMINE  to  estimate  air  entrainment  into  buoyant  balcony  spill  plumes.  It  concerns  the 
numerical  simulation  of  a  reduced  scale  physical  model.  Excellent  agreement  between 


measured  and  predicted  mass  fluxes  is  demonstrated  in  the  regions  of  low  heights  of 
plume  rise,  but  the  two  apparently  diverge  at  greater  heights  of  plume  rise. 

Currently,  the  development  of  a  computational  fire  dynamics  code  is  underway 
at  NIST.  The  code,  Fire  Dynamics  Simulator  (FDS),  including  its  visualization  tool, 
SMOKEVIEW,  is  looked  to  as  the  future  research  tool  for  fire  investigation  |261.  It  is 
based  on  a  Large  Eddy  Simulation  (LES)  with  a  Smagorinsky  turbulence  sub-model. 
The  hydraulics  equations  were  initially  developed  by  Baum  and  Rehm  |271.  The  general 
usage  and  application  of  the  code  has  been  reported  by  McGrattan  et  al.  |2S|.  This  code 
has  been  used  to  study  the  isolated  plume  dynamics  |291.  The  combustion  model  used 
there  is  based  on  simulating  the  combustion  by  releasing  energy  from  Lagrangian 
thermal  element.  The  radiation  is  simulated  through  ray-tracing  from  each  thermal 
element.  A  45  kW  pool  fire  plume  is  simulated.  The  agreement  is  remarkably  good  in 
the  intermittent  flame  and  plume  region,  but  not  in  the  continuous  flame  region.  Since 
the  energy  release  from  the  fuel  particles  was  uniform  in  time,  dense  particles  near  the 
burner  surface  caused  the  temperature  to  be  over-predicted.  This  discrepancy  limits  the 
application  of  the  code  in  predicting  flame  heat  flux,  and  also  in  its  general  use  for  fire 
spread  prediction. 

Based  the  need  for  improvement  in  predicting  flame  heat  flux,  an  alternative 
combustion  model  has  been  realized  in  the  code  by  McGrattan  et  ai.  |3t)|.  The  new 
combustion  model  is  based  on  infinitely  fast  chemistry  and  mixture  fraction.  A  new 
radiative  model  using  Finite  Volume  Method  (FVM)  for  radiation  transfer  was  also 
incorporated.  In  order  to  evaluate  this  new  model  for  general  applications  in  fire 
growth,  it  was  decided  to  investigate  the  axi-symmetric  fire  plume  since  it  was 
extensively  studied.  However,  since  the  experimental  data  are  too  abundant  and  not 
easily  accessible,  a  review  of  existing  correlations  is  made  to  determine  the  most 
representative  and  accurate  approach  to  use  for  these  properties.  A  series  of  simulations 
were  carried  out  in  order  to  determine  these  gross  fire  plume  properties,  e.g.  the 
centerline  time-averaged  Froude  number,  temperature,  and  flame  height,  etc. 

2.  Methodology 

Here  we  briefly  summarize  the  new  models  adopted  in  the  code,  which  is 
supplied  by  McGrattan  et  al. |301. 

2.1.  Hydrodynamic  and  Turbulence  Model 

Details  of  the  hydrodynamics  and  turbulence  model  are  supplied  in  |28]|291.  For 
the  current  study,  the  fundamentals  of  Large  Eddy  Simulation  are  unchanged. 

2.2.  Combustion  Model 

The  new  combustion  model  used  here  is  based  on  mixture-fraction  based 
infinitely  fast  chemistry  kinetics.  The  general  form  of  the  combustion  reaction  is, 

Vo[0]  +  vf[F]-»^vp[P]  (1) 

The  mixture  fraction  Z  is  defined  17711801  as, 
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The  mixture  fraction  satisfies  the  conservation  law 

p——  =  V  ■  pDVZ  (3) 

Dt 

The  assumption  that  the  chemistry  is  "fast"  means  that  the  reactions  that 
consume  fuel  and  oxidizer  occur  so  rapidly  that  the  fuel  and  oxidizer  cannot  co-exist. 
The  flame  sheet  is  the  location  where  fuel  and  oxidizer  vanish  simultaneously. 
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When  dealing  with  a  wide  range  of  fire,  this  definition  is  not  enough  to  produce 
the  correct  flame  height.  Sometimes  the  code  may  not  have  enough  resolution  to  capture 
the  combustion  and  relating  dynamics.  An  effective  Z,  is  proposed  to  help  the  code 

capture  the  combusting  region. 

Zfeff  (  D'\ 

—L^L  =  min  1,C - ,  (4b) 
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Here  C  is  an  empirical  constant  to  be  used  for  all  fire  scenarios  (here  it  is  0.025  in 
FDS2.0),  <5x  is  the  grid  size  and  D*  (=z*  herein)  is  the  local  characteristic  length  of  the 
fire,  used  as  a  reference  length.  As  the  resolution  of  the  calculation  increases,  the 
effective  Z;  e(f  will  approach  the  ideal  value,  Z,  of  equation  (4a).  The  benefit  of  the 

expression  is  that  it  provides  a  quantifiable  measure  of  the  grid  resolution  that  takes  into 
account  not  only  the  size  of  the  grid  cells,  but  also  the  size  of  the  fire. 

The  mixture  fraction  is  a  dimensionless  value  so  it  cannot  decide  the  species 
distribution  alone.  A  state  relation  for  oxygen  is  needed.  The  ideal  state  relation  for 
oxygen  is  introduced  based  on  the  assumption  that  f"?1  and  oxidizer  cannot  co-exist. 

For  Z<Zf,Y0(Z)  =  Y;  (1  -  /%/, )  (5a) 

ForZ  >  Zf ,  Y0(Z)  -  0  (5b) 

The  heat  release  rate  per  unit  volume  is  based  on  Huggett's  relationship  1431  of 
oxygen  consumption,  where  A h0  is  assumed  constant  for  most  fuels. 

q"  =  A  hQm"'  (6) 

where  the  mass  burning  rate  is  based  on  the  oxygen  consumption  rate  given  as, 

-m;-V-(pD^Vz)-^-V-pDVZ  (7) 

The  ability  to  predict  the  VZ  over  the  cell  is  a  reason  for  equation  (4b). 

Specifically  for  our  fire  simulation,  by  default  in  the  code,  we  can  have  the 
following  diagram  showing  the  relationship  between  species  for  fires  of  a  typical  fuel, 
propane.  Her eAh,  =13.1  MW/kg  is  used  for  propane. 
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Figure  I.  Relationship  of  mixture  fraction  and  species  concentration  for  propane 


2.3.  Radiation  Model 

The  Radiative  Transport  Equation  (RTE)  for  a  non-scattering  gray  gas  is 

r  .V/(.v',r)  =K(x)[Ib(x,r)  -I(x,r)]  (8) 

Where  I{x,r)  is  the  radiant  heat  flux,  Ih(x,r)  is  the  blackbody  radiation  intensity,  s  is 
the  unit  normal  direction  vector  and  K(r )  is  the  gray-gas  absorption  coefficient.  The 
radiative  loss  term  in  the  energy  equation  is 

- W-qr(x)  =  /c(x)-[f7(x)-4^/6(jr)],  U(x)  =  I(x,r)dQ  (9) 

The  net  radiant  energy  gained  by  a  grid  cell  is  the  difference  between  that  which  is 
absorbed  and  that  which  is  emitted.  The  source  term  is  defined  as 

f koT  '  l  ji ,  Outside 

kI,  J  (10) 
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Here  q"  is  the  chemical  heat  release  rate  per  unit  volume  and  x,  is  the  local  fraction  of 
the  energy  emitted  as  thermal  radiation.  The  radiant  heat  flux  vector  qr  is  defined 


qr(x)  =  Js/(x,s)dQ  (11) 

Through  the  simulations,  we  use  the  default  parameters  for  radiation  calculation,  since 
we  have  no  better  experimental  result  on  radiation  prediction  in  flames. 


3.  Plume  Theory  and  Benchmark  Correlations 

There  has  been  some  work  done  with  FDS  on  isolated  fire  plumes  .  Our 
purpose  is  to  examine  the  code  for  a  wide  range  of  fire  plumes.  The  former  needs  a 
specific  fire  experiment  while  the  latter  needs  general  experimental  correlations  for 
comparison.  In  the  following,  plume  theory  will  be  used  to  extend  existing  correlations 
to  a  more  complete  and  accurate  representation  of  data.  While  some  specific 


correlations  are  considered,  it  should  not  be  inferred  that  others  not  considered  are 
inferior. 

3. 1 .  Scaling  of  Fire  Phenomena 

The  appropriate  scaling  parameters  for  most  of  the  features  of  fire  plumes  have 
not  been  universally  established  by  all  investigators.  However,  the  scaling  laws  for  the 
size  of  the  diffusion  flames  have  been  developed  by  a  number  of  authors,  and  some 
consensus  has  been  reached  that  a  dimensionless  fire  power  131 
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or  a  dimensional  parameter  that  is  proportional  to  Q’d  ,  such  as  Q!  D5  2  can  be  used  as 
the  scaling  parameter.  Flame-length  data  discussed  by  Zukoski  1181  and  Hasemi  161 
indicate  that  scaling  laws  for  the  visible  flame  length,  and  possibly  for  other  parameters 
of  interest  including  the  entrainment  rate,  have  a  different  functional  form  for  several 
ranges  of  Q*D  ,  and  must  be  considered  separately  for  each  regime. 


It  is  also  pointed  out  by  Heskestad  1321  that  Q*D  does  not  appear  to  be  a  sufficient 
scaling  parameter  for  the  flaming  region.  Flame  height  measurements  data  support  a 
new  scaling  parameter  N. 
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In  practice,  the  distinction  from  Q‘d  is  important  mainly  when  atmospheric  conditions 
differ  significantly  from  standard.  So  here  we  still  rely  on  Q"d  as  the  scaling  parameter. 

According  to  Hasemi  et  al.  |33i,  natural  fires  always  fall  into  the  range  of  Q'd  = 
0.1  to  Ql)=  10.0,  Figure  2  summarizes  the  correlation  between  fuel  size  and  Q\y  for 
typical  natural  fires  in  unconfined  environments.  This  demonstrates  that  the  range  of 
Qd  for  most  building  combustible  objects  with  the  size  of  0.1-2  m  is  approximately 
0, 3-3.0,  whereas  measurements  of  the  plume  properties  have  been  made  mostly  for  a 
wider  range  of  QD .  In  order  to  investigate  the  significance  of  the  flame  geometry  on  the 
plume  properties,  Hasemi  et  al. 161  measured  the  centerline  and  the  radial  distribution  of 
temperature  and  upward  velocity  on  porous  round  burners  with  propane  as  the  fuel  for 
0.4<  Q*d  <4,0.  Here  our  computed  simulation  range  is  0.1  <Q*D  <  10.0,  which  covers  the 
range  of  Q*p  for  most  common  fires.  In  order  to  develop  the  appropriate  theoretical 
basis  for  the  general  dimensionless  correlations,  we  will  examine  specific  fire  plume 
features. 


Figure  2.  Dimensionless  QD  for  most  natural  fires’” 

3.2.  Plume  Froude  Number 

An  inviscid  Bernoulli  equation  analysis  for  the  center  streamline  of  a  buoyant 
plume  with  a  uniform  surrounding  atmosphere  at  T( „  gives, 

+  P  +  Pgz  =  P(0)  (14> 

The  atmosphere  imposes  the  pressure  as, 

p{Q)-p  =  pxgz  (15) 

And  since  the  pressure  is  nearly  constant,  the  perfect  gas  law  gives 

pT  =  (16) 


so  we  have 


This  gives  the  inviscid  steady  value  along  the  centerline  a  buoyant  plume.  In  a 
review  of  real  fire  plume  data,  Beyler  1121  summarized  many  experimental  results  and 
found  the  plume  centerline  Froude  number  ranged  from  about  1.16  to  2.0  for  a  wide 
range  of  data  with  an  average  of  about  1.5  in  the  non-combusting  and  flame  intermittent 
region.  Quintiere  and  Grove  1171  have  shown  that  F/-  =  l,5±0.1  for  both  line  and  axi- 
symmetric  fire  plumes  along  the  centerline  in  their  entire  domain,  combusting  and  non¬ 
combusting  regions.  Heskestad  1321  uses  an  alternative  parameter  defined  as 
£=(9.lf  -Fr'12  in  the  non-combusting  plume  region.  This  value  is  found  to  be  2.2 


(Heskestad  |321),  2.7  (McCaffrey  |41),  2.4  (Kung  et  al.  |tD|)  corresponding  to  a  range  of  Fr 
from  1.29  to  1.94. 

It  appears  that  a  multitude  of  data  from  fire  plumes  suggest  a  fairly  constant 
time-averaged  Fr  along  the  centerline.  Since  it  ties  the  velocity  directly  to  the 
temperature  change,  it  is  a  valuable  feature  and  useful  tool  in  fire  induced  flow  research. 
The  approximate  constancy  of  Fr  also  suggests  a  small  and  uniform  departure  of  real 
fires  from  inviscid  flow. 


3.3.  Centerline  Temperature  and  Velocity 

Dimensional  analysis  of  the  governing  equations  can  show  that  the  centerline 
plume  temperature  satisfies 

f(Xr,z/z\D/z *)  (18) 


where 


Q 


P«C 


,J~y[g 


(19) 


is  a  characteristic  plume  length  scale,  and  D  is  the  effective  diameter  of  the  fuel  source 
(the  diameter  of  the  circle  with  an  equivalent  area).  Theoretical  analysis  1161  can  show 
the  general  power-law  form  of  this  relationship  in  the  combusting  and  the  non- 
combusting  regions  of  the  plume.  The  coefficients  in  these  relationships  must  be  set 
from  experimental  data. 

In  the  non-combusting  plume  region,  Heskestad  141  developed  a  relationship  in 
terms  of  a  virtual  origin  zn,  offset  below  the  base  of  the  fire,  to  account  for  the 
dependence  on  D: 
i.e.  for  z  >  z  f  , 


T-7L 


=  CT{  \-Xf) 


13  Z 


-5 13 


(20) 


where 


^  D 

—  =  -1.02— +  1.37  (21) 

z  z 

gives  the  virtual  origin.  Heskestad  1141  obtained  CT  =9.1  for  the  domain  of  pool  fires, 
0.04  2.0, =  (;*/Z))V".  Quintiere  and  Grove  1171  found  CT  =  10.6  based  on  far 

field  data  of  Yokoi |21,  and  this  variation  suggests  the  current  level  of  uncertainty  in  CT  . 

From  the  constancy  of  Fr,  the  centerline  dimensionless  velocity  can  be 
consistently  written  as 


r-C,( l-xj 


111  2  ~  Z(t 


(22) 


where  Cv  is  3.4  according  to  Heskestad  "4|,  or  4. 17  from  Quintiere  and  Grove1'71 


For  the  combusting  region,  it  is  assumed  that  the  entrainment  rate  of  air  into  the 
plume  controls  the  burning  rate.  With  the  similar  turbulent  mixing  effects,  the  fire 
consumes  the  same  fraction  of  stoichiometric  air.  The  entrained  air  is  about  10  times  of 
the  stoichiometric  air  |l71.  A  theoretical  analysis  based  on  these  assumptions  leads  to  a 
constant  centerline  temperature  in  the  fire  region.  Accordingly, 

/  (23) 
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where  CTf  is  a  constant  dependent  on  a  fit  to  data.  After  checking  Table  1  for  selected 
centerline  flame  temperatures,  chosen  for  accuracy  and  fuel  range,  an  estimate  of  CT!= 
0.50  is  a  good  choice.  For  the  propane  fuel  in  our  simulation,  the  temperature  is 
estimated  as  1041°C  in  the  flaming  region. 


Table  1.  Maximum  accurately  measured  centerline  flame  temperatures. 


Fuel 

Mc  lsa 

kJ/g  air 

D 

(m) 

Q 

(kW) 

Source 

r 

TJ 

1260 

Nat.  Gas 

2.91 

0.3 

17.9 

0.15 

Smith  et.  al. ,33' 

0.50 

1100 

Nat.  Gas 

2.91 

gif 

70.1 

0.25 

Smith  et.  al. 1331 

0.48 

930 

2.96 

wsm 

7713 

0.43 

0.52 

990 

Methanol 

3.07 

1.74 

973 

0.18 

0.59 

This  result  should  serve  as  an  approximate  benchmark  in  the  combusting  region 
to  assess  the  accuracy  of  both  predictive  and  experimental  results.  The  radiation  fraction 
is  a  significant  factor,  and  it  probably  represents  the  most  complex  phenomena.  The 
value  of  X r  is  difficult  to  measure  and  will  vary,  depending  on  fuel  type  and  fire 
diameter.  A  fundamental  predictive  model  for  X r  would  require  more  detailed  data  and 
be  capable  of  resolving  distributions  for  flame  properties. 


3.4.  Flame  Height 

Numerous  correlations  exist  for  flame  height,  and  some  are  summarized  in  the 
SFPE  Handbook  |l31.  Usually,  flame  height  correlations  are  produced  under  a  specific 
situation  and  offer  independent  results  from  each  other.  Among  them,  Heskestad's  |S| 
expression  is  easy  to  use  and  Zukoski's  1261  expression  is  based  on  consistent  data. 
Heskestad  gives: 

z  .  O*2'5 

—  =  15.6- — — - 1.02  (24) 

D  (bhc/sa) 


which  is  commonly  presented  as 

~z~  =  3-88  •  Ql2li  -  1 .02 ,  for  Ah  /  =  2.975  kJ/g  . 

D  L  “ 


(25) 


Based  on  the  three  data  regimes  characterized  by  Qo  ,  Zukoski  developed 


f/_40.o*2,  for  0.03*£;<0.15 

D 

(26a) 

-  3.3-Qq1'3 ,  for  0.15  1.0 

(26b) 

—  =  3.3  •  C?/>2  5  >  for  1.0  <Ql  <40.0 

D 

(26c) 

Both  correlations  (equation  (24)(26))  ignore  the  effect  of  flame  radiation. 
Hamins  et  al.  1111  attempted  to  include  this  effect  with  an  empirical  inclusion  of  Xr 
They  obtained 


—  =  18.5 
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This  expression  does  not  contain  the  correct  dependence  for  X r  according  to  the 
theory  of  Quintiere  and  Grove  |[7).  They  theoretically  derived  the  effect  of  radiation  and 
found  an  implicit  expression  for  flame  height. 

Qd  _ _  (28) 


(X) 

1/2 
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For  z ,  ID  >  0.1,  the  LHS  can  be  approximated  as 


(l-Xj12 

[  0.775  +  0.463—1 

(  D) 


■ ,  giving 


—  « 16.8 
D 


e;2 


15 


(A  hji 


3/5 


-1.67 


(29) 


cpTx 


(i  -x,r 


This  has  the  form  of  Equations  (24)  and  (26).  It  would  appear  that  radiation  does  not 
have  a  strong  effect  on  flame  height.  However,  it  has  the  opposite  trend  in  Hamins’ 
result.  Qualitatively,  as  X r  increases,  the  flame  temperature,  velocity  and  entrainment 
rate  should  decrease,  and  therefore  the  flame  height  should  increase  as  more  air  is 
available  for  combustion. 

For  z,  ID  <  0.1,  the  LHS  of  Equation  (29)  can  be  approximated  as  (z f  ! D)V1  , 
therefore 

—  *■»  2.873x  104 - % - 

D  Uhls] 


\  f  /’^° 


0-^) 


or  using  Ahc  I  sa  =  2.975  kJ/g 


(31) 


Sr- 
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*2 

1-X^’ 


with  the  same  form  as  Zukoski  (equation  26c),  but  now  including  X 


The  range  of  flame  height  data  is  about  0.05  <.  Q'D  <.  100.  All  of  these  formulas 
give  good  accuracy  above  Q*D=  1.0.  There  is  a  larger  disagreement  among  the  formulas 
and  a  larger  scatter  of  the  data  below  Q*D  = 1.0.  We  shall  use  Equations.  (24),  (26),  and 
(29)  as  good  representations  of  flame  height  (shown  in  Figure  3).  In  practice,  we  usually 
take  the  following  parameters:  Ahc!sa  =2.975  kJ/g,  cp  =1.0  kJjkg-K  ,TX  =298  K. 


10 1  10°  10' 

Dimensionless  fire  power,  QD 

Figure  3.  Comparison  of  Flame  Height  Correlations. 


3.5.  Entrainment  number  at  flame  tip 

There  is  an  important  constant  in  relating  the  entrainment  rate  and  flame  height 
prediction.  It  is  an  empirical  number,  which  has  different  values  among  different 
researchers  due  to  different  test  methods  or  analyzing  theories  applied.  This  number,  n, 
is  the  ratio  of  entrained  air  up  to  the  flame  tip  with  the  stoichiometric  air  needed. 

Some  early  research  work  on  fire  plume  was  performed  by  Yokoi  '2|.  Roughly  12 
times  stoichiometric  air  at  flame  tip  was  deduced  by  McCaffrey  and  Cox  1451  from  his 
temperature  and  velocity  profiles. 

In  order  to  fit  his  combustion  plume  model  to  a  large  quantity  of  flame  height 
data,  Steward  1461  found  that  only  400%  excess  air  was  entrained  up  to  the  flame  tip,  i.e., 
the  flame  has  entrained  about  5  times  the  stoichiometric  requirements. 

Using  k-e  -g  model,  Tamanini  |1X|  produced  9  times  stoichiometric  at  flame 
tip.  This  value  was  thought  to  be  low,  due  to  the  deficiencies  in  the  treatment  of  large 
scale  eddies  near  the  point  of  transition  to  turbulence. 

Zukoski's  151  hood  collection  technique  produced  an  equivalence  ratio  of 
0.067  ±  10%  measured  value  at  flame  tip,  which  yields  about  15  times  stoichiometric. 


Based  on  point  measurement  techniques  coupled  to  an  integral  model  of  the 
fluid  dynamics  behavior,  McCaffrey  and  Cox  1451  got  the  number  to  be  about  20  times 
the  stoichiometric  requirement. 

me  f  =  0.0070  «  20 (sam  f )  (32) 

Beyler  |4/(  measured  the  entrainment  rate  using  the  hood  technique.  He  got 

me  J  /  Q  =  0.045  kg/k/kW,  (33) 

representing  about  13  times  the  stoichiometric  air  requirement  of  the  flame. 

Based  on  Zukoski's  data  on  entrainment  at  flame  tip,  Heskestad  1481  assumed 
self-preserving  density  deficiency  profiles  instead  of  the  assumption  of  self-preserving 
excess  temperature  profiles  used  by  Cetegen  et  al. 1491  It  is  found  that  12  times  the  mass 
stoichiometric  air  is  required  at  the  flame  tip. 

</  =0.00540,  =0.0038(samf)(Ahc/sa)~\2(samf)  (34) 

Using  an  apparatus  similar  to  Ricou  and  Spalding  |501,  Delichatsios  and  Orloff 1511 
found  that  the  air  entrainment  is  about  10  times  the  stoichiometric  requirements  at  the 
flame  height  of  turbulent  buoyant  fires.  This  number  had  been  used  directly  to  derive 
the  flame  height  correlations  from  the  entrainment  correlations  |9). 

Examining  a  lot  of  data  using  an  integral  model  approach,  Quintiere  et  al.  1171 
derived  the  entrainment  and  the  flame  height  separately,  and  by  comparison  of 
coefficients,  found  n=9.6  for  both  axi-symmetric  and  linear  plumes. 

Hence,  the  range  of  n  is  found  to  be  10±  5 . 

4.  Simulation  Results 

4.1.  Simulation  Summary 

Generally  in  fire  research,  a  circular  burner  is  used  more  widely  than  a  square 
burner  which  is  used  in  McCaffrey's  experiments.  Measurements  by  Hasemi  et  al. ,6f  for 
square  and  round  gas  burners  show  little  effect  of  the  burner  shape  when  for  the  square 
burners,  D  is  taken  as  the  diameter  of  a  circular  burner  of  equal  area.  Subsequent 
measurements  of  Hasemi  and  Nishita  1351  have  indicated  tendencies  to  slightly  higher 
flame-height  ratio  than  circular  burner  data  for  small  fires.  Investigation  by  Orloff  13611371 
revealed  that  the  height  of  a  peripheral  rim  or  lip  above  the  surface  of  a  circular  pool 
fire  affects  the  burning  rate  and  flame  shape,  but  has  little  effect  on  the  flame  height 
(except  through  Q‘d).  Other  geometric  effects  on  flames  have  attracted  attention. 
Without  further  comment,  these  include  interactions  with  ceilings,  interaction  with 
walls  and  corners,  and  interaction/merging  of  proximate  flames. 

Here  for  the  numerical  plume  simulation,  we  are  dealing  with  an  unconfined, 
free  burning,  square  burner  surface,  no  rim  or  lip,  pool  fire  with  propane  as  typical  fuel, 
and  with  fire  sizes  ranging  from  Q*D  =0.1  to  10,0,  (For  our  0,3x0. 3m  square  burner,  the 
fire  size  ranges  from  7.37  kW  to  737  kW),  Unless  stated,  all  simulations  are  based  on 
constant  radiative  loss  fraction  ( X r  =0.35)-  The  mass  flow  rate  of  the  propane  was 

specified  and  supplied  with  uniform  velocity.  The  platform  for  the  simulation  is  the 


Pentium  IV,  2.0GHz,  IGB  RAM,  and  Windows  2000.  All  cases  are  run  for  10  seconds 
of  simulation  with  last  4  seconds  for  steady  state  averaging. 


Table  2. 


ist  of  simulation  inputs 


Q* 

Grids 

0(  input) 

7* 

Grid 

R* 

Q(out) 

Error_Q 

Run_time 

kW 

m 

cm 

kW 

hrs 

0.10 

64x64x64 

7.37 

0.13 

0.69 

0.0510 

7.33 

-0.01 

3.91 

0.20 

64x64x48 

14.74 

0.18 

0.91 

0.0510 

14.49 

-0.02 

2.64 

0.50 

48x48x64 

36.86 

0.26 

1.25 

0.0487 

37.41 

0.01 

2.29 

0.80 

32x32x96 

58.98 

0.31 

1.56 

0.0505 

58.25 

-0.01 

1.39 

1.00 

32x32x80 

73.72 

0.34 

1.69 

0.0499 

74.84 

0.02 

1.23 

1.50 

32x32x80 

1 10.59 

0.40 

2.00 

0.0502 

114.79 

0.04 

1.19 

2.00 

32x32x96 

147.45 

Hi 

2.25 

0.0504 

0.06 

1.30 

5.00 

32x32x96 

368.62 

3.25 

0.0504 

0.12 

1.08 

8.00 

32x32x96 

589.79 

mmSm 

3.88 

0.0498 

652.86 

0.11 

1.02 

10.00 

32x32x96 

0.85 

4.25 

803.10 

0.09 

0.30 

48x48x144 

fmmrn 

0.21 

1.05 

21.54 

As  Zukosift1  pointed  out,  when  Q'd  <0.10,  the  combustion  region  is  divided 
into  many  flamelets  that  lean  toward  the  axis  of  the  source  and  that  have  flame  lengths 
much  smaller  than  the  diameter  of  the  source.  For  values  above  1 .0,  the  flames  are  more 
than  three  diameters  high  and  have  a  columnar  form.  In  the  middle  regime  the  short 
flamelets  merge  as  QD  increases  and  form  a  single  short  flame  near  QD  =0. 1 .  Here  the 
simulations  produced  similar  results  as  shown  in  Figure  4. 


Figure  4.  Snapshot  of  fire  plume  at  QD  =0.20  and  QD=\Q.O 


The  reason  that  we  could  not  show  results  for  Q*D  <  0.2  fires  here  is  that  for  very 
low  fire  intensity,  the  Z  =  Z,,  contour  will  be  very  close  to  the  burner,  even  for  a  very 
well  resolved  gridsize.  As  the  cells  get  larger,  the  mixture  faction  value  at  the  burner 


surface  is  too  small  to  reach  its  stoichiometric  value.  Sometimes,  for  low  enough  Q*D  , 
the  Z  =  Z  will  be  located  inside  the  burner,  which  is  invisible  in  SMOKEVIEW1521, 

the  visualization  tool  for  FDS.  We  still  can  see  the  flame  structure  if  a  smaller  mixture 
fraction  value  is  used,  but  this  has  less  physical  meaning. 

4.2.  Flame  Height 

There  are  some  different  definitions  of  flame  height.  The  averaged  luminous 
flame  height  is  generally  defined  as  the  height  at  which  the  flame  is  observed  at  or 
above  that  height  50%  of  the  time  |l61.  Hasemi’s  flame  height  was  defined  as  the  limit 
above  which  flame  tips  have  never  been  observed  by  eye  |6L  Visual  observation  tends  to 
yield  slightly  over-estimates  of  the  50%  flame  heights  |l6].  Video  image  analysis  seems 
to  provide  the  best  estimates  though  averaging  a  number  of  one-second  exposure 
photographs.  Cox  and  Chitty  found  the  flame  lengths  from  3  seconds  time  exposure  by 
photographs  of  fires  171 . 

For  our  simulation,  there  are  many  parameters  available  that  can  serve  as 
indicators  of  flame  height,  such  as  temperature,  fuel  concentration,  mixture-fraction,  or 
energy  release,  etc.  Here  we  couldn't  use  50%  intermittency,  as  the  frequency  of  fire 
pulsation  is  still  not  accurately  resolved.  Instead,  we  choose  a  more  fundamental  one, 
local  energy  release  rate  (Q'"),  as  the  parameter  for  flame  height.  For  calculating  the 
flame  height,  first  we  average  the  recorded  slice  file  for  Qm  .  Generally,  we  choose  to 
average  the  frames  of  4  s  steady  fire  simulation.  Then  the  Q  "  value  was  integrated  over 
the  combustion  domain  upwards  from  the  burner  surface.  So  we  have  an  energy  release 
rate  as  a  function  of  height.  The  critical  value  for  flame  tip  was  considered  as  95%, 
99%,  or  99.9%.  After  a  series  of  simulations,  we  decided  that  99%  is  a  good  indicator 
for  flame  height  (Figure  5.).  That  is  to  say,  our  flame  tip  is  defined  as  the  level  below 
which  99%  of  the  chemical  reaction  occurred. 

The  impact  of  gridsize  on  simulation  results  is  an  important  consideration  for 
any  CFD  analysis.  Generally,  good  CFD  software  should  have  a  relative  constant  result 
in  a  certain  gridsize  range.  Baum  et  al.  1291  suggested  that  to  correctly  capture  the  impact 
of  burner  surface,  it  is  necessary  to  have  10x10  grids  over  the  burner  surface.  Here  since 
we  are  dealing  with  a  wide  range  of  fire  intensity,  this  empirical  rule  is  not  enough.  The 
choice  of  gridsize  should  reflect  not  only  the  impact  of  fire  dimension,  but  also  the  fire 
size.  The  characteristic  length,  defined  by  equation  (19),  combines  the  impact  from  both 
burner  diameter  and  the  fire  size,  so  it  is  widely  used  in  fire  research  field.  Here  we  can 
define  the  resolution  of  fire  plume  simulation  as  a  dimensionless  parameter, 

R-. *  _  max(&,3y,&)  (35) 

z  * 

To  simulate  the  fire  correctly,  non-uniform  grid  is  allowed  and  the  aspect  ratio 
of  grid  size  should  not  exceed  2.  Here  we  have  uniform  grid  size  in  all  three  dimensions 
in  all  runs,  i.e.  5x  =3y  =Sz  .  The  simulation  resolutions  for  all  cases  are  listed  in  Table 
2. 
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Figure  5.  Flame  height  prediction  as  a  function  of  combustion  termination. 

After  a  lot  of  simulation  runs,  the  optimum  resolution  for  flame  height 
simulation  was  found  to  be  R*  =0.05.  Above  this  value,  the  simulations  tend  to  under¬ 
predict  the  flame  height.  Below,  then  the  flame  height  is  over-predicted.  It  is  found  291 
that  the  plume  dynamics  can  only  he  accurately  simulated  if  the  resolution  limit  is  about 
R*=0. 1  or  smaller.  Our  result  is  even  strict.  The  flaming  region  has  the  chemical 
reactions  and  wide  range  of  turbulence  scales,  so  a  finer  grid  size  is  needed. 

A  closer  look  at  other  critical  parameters  at  the  flame  height  (calculated  using 
Zukoski’s  correlation)  shows  the  following  result  (Figure  6.). For  fires  between  Q'd  =  1.0 
to  £^  =  10.0,  the  temperature  at  flame  height  is  about  T=350  “C,  which  is  close  to 
Zukoski’s  value.  The  critical  mixture  fraction  Z=0,01  is  also  close  to  the  result 
calculated  using  equation  4(a).  For  small  fires  (£9^,  <0.5)  the  reaction  happens  only  one 
or  two  grid  cells  from  the  burner  surface,  the  high  value  at  flame  tip  is  caused  mainly  by 
diffusion  alone,  where  the  turbulence  is  low  and  the  entrainment  is  not  enough.  The  fire 
is  more  laminar-like,  rather  than  turbulent.  But  the  temperature  prediction  follows  the 
correlation  over  a  wide  range  of  fires. 


V  /  A 


Figure  6.  Temperature  and  Mixture  fraction  at  the  flame  tip. 
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As  Tamanini  1441  pointed  out,  since  our  combustion  model  is  based  on  the 
concept  of  mixture  fraction,  the  temperature  and  the  mixture  fraction  share  the  same 
governing  equation.  When  using  a  constant  radiative  loss  fraction,  the  temperature  and 
the  mixture  fraction  are  linearly  related,  so  these  two  lines  should  converge  into  one 
line,  or  parallel  to  each  other.  Figure  6  shows  this  trend.  As  mixture  fraction  level  is  an 
indication  of  the  mixing  level,  it  is  also  deduced  that  the  mixture  fraction  level  at  flame 
tip  is  0.013,  which  is  about  one  sixth  of  its  stoichiometric  value.  This  shows  that  the 
flame  entrains  about  6  times  stoichiometric  air  to  the  flame  tip. 

Currently  there  are  no  direct  data  available  on  temperature  measurement  at  the 
flame  tip.  But  using  McCaffrey's  well-accepted  temperature  distribution  |l3<  and 
Heskestad's  simple  flame  height  correlation  |H1,  we  can  deduce  an  empirical  value  of 
temperature  at  the  flame  tip.  McCaffrey's  correlation  for  temperature  has  the  form, 
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(36) 


Here  the  flame  tip  should  be  somewhere  in  the  intermittent  region.  The  according 
constants  for  intermittent  region  are,  k=1.9,  C=0.9,  and  n=0,  The  flame  height  using 
Heskestad's  correlation,  which  is  very  simple  in  form: 


zf  =0.235g-1.020 


(37) 


The  solid  line  in  Figure  6  shows  the  temperature  prediction  from  equation  (36)  and  (37). 


4.3.  Buoyant  Froude  Number 

Figure  7  shows  the  computational  results  for  the  averaged  centerline  Froude 
number.  The  Froude  number,  after  the  flow  develops  away  from  the  burner  surface, 
ranges  from  2  for  a  very  small  fire  to  1.1  for  a  very  large  fire.  For  fire  intensity  Q*D 
from  0.2  to  1.0,  which  covers  McCaffrey  141  experimental  range,  the  Froude  number  lies 
within  the  range  of  1.5  to  1.8,  while  McCaffrey  got  1.62.  This  value  is  virtually 


invariant  throughout  the  fire  plume  region,  from  the  intermittent  zone  right  down 
through  the  continuous  flames  at  the  base  of  the  burner. 


For  low  Qd  fires,  the  reaction  zone  is  very  small,  the  flow  is  less  turbulent  and 
perhaps  inviscid  and,  i.e.,  Froude  number  is  close  to  2.0.  For  fires  of  large  Q*D  ,  the 
flame  is  more  jet-like.  The  computed  results  are  in  the  range  of  experimental  results,  i.e. 
Fr  =  1.5  ±  0.1,  and  show  a  systymmatic  effect  of  Q*D  which  appears  related  to  viscous 
effects. 

4. 4.  McCaffrey 's  Experiment 

McCaffrey  |4i  did  a  series  of  classical  experiments  in  1979.  He  used  natural  gas 
at  varied  controlled  rates  on  a  0.3  m  square  burner.  His  experimental  data  is  used  widely 
to  test  combustion  models.  Here  we  directly  compare  the  simulation  result  with  his 
experiments.  His  Q'd  is  0.3  for  this  fire. 


Figure  8.  Temperature  simulation  (left)  with  McCaffrey’s  22  kW  fine  (right)14' 

From  Figure  8,  it  is  apparent  that  the  vertical  temperature  distribution  is  fairly 
good,  while  the  flame  is  slightly  thinner  than  the  real  measurement.  A  caution  should  be 
made  here  that  the  temperature  in  McCaffrey's  experiment  is  acquired  without 
correction  of  thermocouple  data.  So  the  numerical  simulation  is  still  low  for  this  fire 
size  (approx.  Q*0=  0,3).  The  centerline  temperature  distribution  shows  this  trend. 


Figure  9.  Centerline  temperature  distribution  with  McCaffrey's  22  kW  fire |4) 


Deep  in  flames,  the  mean  temperature  reaches  approximately  900  "C.  Values  of 
mean  temp  as  high  as  1000  "C  have  been  observed  by  Gregory  1401  in  large  9x1 8m2  pool 
fires.  Also  Smith  et  al  1,11  recorded  a  peak  mean  temperature  of  1260  "C  for  17.9  kW 
natural  gas  flame.  The  simulated  temperature  fits  well  in  the  plume  region,  while  over¬ 
predicted  in  the  flaming  region.  The  velocity  from  McCaffrey's  experiment  contracted 
together  into  one  curve,  our  simulation  fits  well  with  the  measurement  data. 

4.5.  Temperature  Distribution 

Figure  1 1  shows  the  temperature  distribution  along  the  centerline  for  a  range  of 
fire  size.  It  shows  that  the  simulated  temperature  distribution  follows  the  correlation 
well  in  plume  region.  But  the  temperature  rise  in  the  combusting  region  is  still  a  little 
high,  and  the  temperature  drops  earlier  than  expected  in  the  plume  region.  For  the 
temperature  near  the  burner  surface,  the  unexpected  low  values  for  large  fires  and  the 
high  values  for  small  fires  come  from  the  main  assumption  that  the  combustion  is 
assumed  to  be  infinitely  fast  reaction,  so  most  reaction  happens  close  to  a  flame  sheet. 
The  location  of  the  flame  sheet  determines  the  flame  structure  and  all  other  global 
parameters. 


Figure  1 1.  Centerline  temperature  distribution  with  Heskestad's  correlation, 
equation  (20)(21)(23). 

We  do  not  have  good  temperature  data  in  the  flame  over  a  wide  range  of  Q*D  . 
The  correlation  of  Eq.  (23)  is  approximate.  Flere  from  Table  2  and  Figure  5,  it  can  be 
shown  that  there  are  more  grid  cells  over  the  burner  diameter  for  Q*D -0.5,  while  the 
same  number  of  grid  cells  is  over  the  flame  height.  It  is  not  clear  why  the  simulation 
gives  such  a  unusual  effect  in  the  combustion  region,  but  good  results  in  the  far-plume 
region.  It  is  due  in  part  to  a  high  fuel  (cool)  flow  rate  for  high  Q*D  that  lowers  the  flame 
temperature  near  the  burner. 

4.6.  Entrainment  Number 

As  pointed  out  by  Delichatsios  '9|,  the  entrainment  number  n,  the  entrained  air  to 
stoichiometric  air  required  ratio  at  the  flame  tip,  is  important  in  determining  the  flame 
height.  Due  to  the  complexity  of  entrainment  measurements  and  its  variations  in 
theoretical  development,  this  number  is  found  to  have  many  values.  Here  we  check  the 
entrainment  behavior  of  the  flame  height  cases  and  find  the  computed  entrainment 
behavior  at  different  levels  (shown  in  Figure  12).  The  entrainment  numbers  at  flame  tip 
are  shown  in  Figure  13,  where  the  stoichiometric  air  to  fuel  mass  ratio  for  propane  is 
found  to  be  17.3. 
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Flame  spread  and  burning  rate  is  examined  in  a  5-m  high  PMMA  vertical  slab, 
using  the  Fire  Dynamics  Simulator  (FDS),  a  Computational  Fluid  Dynamics  model. 
Coupled  descriptions  of  the  controlling  mechanisms  of  hydrodynamics  model  for 
buoyancy-driven  flow,  mixture-fraction  model  for  combustion,  convection  model, 
banded  gray-gases  transportation  model  for  radiation,  1 -dimensional  conduction  model 
for  pyrolysis  and  burning  rate  model  are  included.  The  approaches  for  pyrolysis  and 
burning  rate  model  are  evaluated  in  a  sub-model  level,  and  the  refined  algorithms  are 
presented. 

The  models  are  first  tested  on  a  non-spread  case  to  examine  the  accuracy  of 
prediction  for  heat  flux  and  gas  temperature.  Results  were  validated  with  the  literature 
and  the  radiation  sub-model  is  tested  with  various  parameters. 


Further  predictions  are  made  on  a  more  complex  large-scale  flame  spread  case 
involving  a  PMMA  wall,  5  m  tall  and  0.6  m  wide,  with  perpendicular  steel  draft 
containing  the  flame  sheets.  Predictions  on  heat  flux,  gas  and  solid  temperature,  flame 
spread,  and  mass  burning  rate  are  evaluated. 

Generally,  good  agreement  for  flame  spread  and  energy  release  rate  is  achieved 
with  experiments  at  levels  of  spatial  resolution  that  are  compatible  with  the 
computational  demands  of  field  model  predictions  in  compartment  fire  scenarios. 


ACKNOWLEDGEMENTS 


Firstly,  I  would  like  to  thank  my  advisor,  Professor  Jim  Quintiere,  for  his 
continuous  motivation,  patient  guidance  and  sparkling  ideas.  I  am  grateful  to  Dr.  Fred 
Mowrer  and  Dr.  Andre  Marshall  for  serving  in  my  review  committee.  A  special  thanks 
is  extended  to  the  faculty  and  staff  of  the  Department  of  Fire  Protection  Engineering  at 
the  University  of  Maryland.  Specifically,  I  would  like  to  thank  Dr.  Milke,  Dr.  di  Marzo 
and  Dr.  Trouve  for  instruction  along  the  way.  I  cannot  thank  Cindy  Gilbert,  Allison 
Spurrier  and  Sylvia  Keane  enough  for  the  help  and  support  provided.  There  are  many 
others  whose  friendship  and  support  have  been  important  to  me  during  the  past  2  years. 
A  partial  list  follows:  Tingguang,  Alwin,  Boonmee,  Kauro,  Ali,  Clarence,  Johnny, 
Yao.. .  Thank  you  very  much!  It  is  truly  an  honor  to  have  received  an  education  from 
the  Department. 

This  project  would  not  have  been  possible  without  the  financial  support  of  the 
National  Institute  of  Standards  and  Technology  monitored  by  Dr.  Thomas  Ohlemiller. 
In  particular,  I  would  like  to  thank  Dr.  Thomas  Ohlemiller  and  Dr.  Kevin  McGrattan  of 
NIST  for  their  infinite  patience  and  invaluable  advice  while  working  on  this  project. 

Finally,  words  alone  cannot  express  the  thanks  I  owe  to  my  wife,  Helen,  for  her 
endless  love,  encouragement  and  assistance. 


1 


TABLE  OF  CONTENTS 


Ll SI Ol  i  Kill l<]  S . V 

Li  SI  ( )l;  [ADLLS . vii 

NOMENCLATURE . viii 

1  INTRODUCTION . 1 

1  .  I  Ml'l  l\  ,H  )Oll  1 

1.2  Background  6 

1 .  '  Approaches  and  Objectives  8 

2  1  1 K I  ■  1)Y  N.W1K  S  SI  Ml  I  A  I  OR . 10 

2.1  UynyTal  andj  Ivdrodynamics  Model  10 

2.2  l)iiHisi\e  I  e  in  is  <1.1,  S  j  11 

2J  Combustion  Models  13 

2.4  (  onvecln  e  I  leal  I  lanslci  to  W  alls  18 

2.5  Radial; ye  I  leal  I  rnnxlcr  to  W  alls  18 

2 .  ('  I ’y  rolys  is  Model  Thermally  i  lock  Solid  21 

2  A  Proposed  Algorithm  lor  Surlacc  I  cniperatnre  ol'  Thermal  I  v-lhick  Material  23 

2.N  Proposed  Algorithm  lor  Burning  Rate  ol  1  herniallv  I  hick  Material  24 

3  BA(  K  S  U  All  l  IR1  MSI  S . 30 

5.1  I  )cscnption  ol  the  W  all  !  ire  1  \pci  imeni  30 

3.2  Back’s  1  xpei  imcnlal  Results  31 

4  I  A  All  A  I  ION  Oi  I  III  l!)S:  PR  I  DICTION  ON  HACK’S  \\  A I  1  TEST . 34 


it 


4.1 


IDs  Input 


34 


35 


1.2  I  DS  Predictions  on  t he  (  enterhne  I  enipcrature 

4.3  I  DS  Predictions  on  (lie  W ail  I  leal  I  lux  I  )ist ?  ib.it u>n  36 

4.4  IDS  Predictions  on  the  Haute  llci.ieht  38 

4.5  Summary  ot  the  I  DS  Predictions  on  Heck's  VVnll  hire  I  csls  39 

5  1  M' S  |,AR(il-.  St_Al  1'.  I’M  VI A  W  A l.l  1  IK  1  I  I  ST . 40 

5.1  I  V1R(  5_0l)-kW- Settle  Llammahihty  Apparatus  ami  the  1  ,\penmcnlal  Setup  40 

5.3  AC  lose  Look  at  I'MYSAA  1  hernial  Proper!  ics  44 

5.3  I  VI RC  I  - . x per i menta I  Results  49 

P  1  \  Ml  A  l  ION  Oh  I  111  Ills  PR  |  .DlC  I  ION  ( >\  II II  PM  VIA  \\  A 1 .1  1  I  S  1  ...51 

L  I  I  I  )S  Implementation  lssik>  51 

Ci. 3  C  omparison  ol  the  Ileal  Release  Rate  I  hstoi y  54 

()..3  e  oni  pan  si  »n  ol  the  1  lame  Spread  55 

0.4  (  omparison  ol  the  lotal  l  leat  I  lu\  57 

0.5  C  pniptinson  ol  the  Hut  nine.  Rate  58 

7  i  VALDA  I  I()NJ)I  III]  IDS  PR  h  Die  I  ION  ON  I  III.  SV1AI.1  SCALP.  I’M  VIA 
WALL  I  1ST . 61 

7. 1  Small-scale  I  lame  Spread  I  xpenmemtal  Setup  61 

7.3  SOI  II'  .ind  LDS  ( Calculations  61 

7  J  I  he  PM  VI  A  I  hernia i  Properties  62 

7.4  C  omparispn  ol  the  Results  65 

S  SUM VI A R Y  AN IHONLI.  I : SION'S . 71 

')  RL  ILRINC  LS . 77 


ill 


iv 


LIST  OF  FIGURES 


l  injure  I  R  elaUonship  ol  mi  x  lure  I  hi  chon  aiul  species  concentration  lor  burning 


PMMA . 16 

Figure  3  J_.lv>11  Mass  I  pansier  Processes  lor  a  I  hermajjv  I  hick  Material . 26 

Figure  3  Experimental  Setup  loi  Meek's  Wall.  I  nc . 31 

Finn  re  3  Measured  centerline  heal  Hun . 33 

l  igurej'  l-'l )S  Simulation  I ).oniain  lor  Mack's  Wall  Fire  Test . 34 

Figure />  (  cnierlmc  ( ins  1 cnipcralure  (1 .25  cm  against  wall..}  \  1 1  !  esls . 36 

Figure.7  _  I  lent  _F]ux  Disirihutioii . 37 

Figure  x  Ileal  Flux  Distribution  lor  l  est  (> . 38 

Figure  l>  Maine  I  height . 39 

Figure  111  ._  Full  Scale  i'MMA  Wall  I  ire  lest . 4] 

Figure  1  1  I  MR(  5(>-kW  I  lanimahiliiv  Apparatus  as  shown.  I  lie  5()()-kW  apparatus  is 

sun  i  I  at  except  lot  s  i/e . 42 

Figure  12  PMMA . 45 

Figure  13  Heat  l  apacitx  ol  \  artous  sol  id.  and  molten  I’MMA . 46 

Fjgure  14  I  hernial  C  onductiviK  ol  \  at  ions  Solid  and  Molten  PMMA . 47 


jgure  15  I  hermoeouple  Reading  on  the  Surface  ol  the  Samples  at  the  onset  of 


pyrolysis . 48 

Figure  16  i  otal  Heat  I  lux  I  tistnhui ton  at  Vat  tons  j  line .  49 

Figure  I7  Pyrolysis  Height  History  . 50 

Figure  IN  I  leat  Release  Rate  I  list  pry . 50 


\ 


Heme  l1-*  Inpul  f  ile  lor  i  I ) S  Simulation  ol  the  I  VI  Wall  I  ire  i  esi . 53 

I ■  inure  2 < I  (,  oi  npansi > 1 1  ol  1 1 R  R  U  i si ore  .  A  Iter  Syn cl i rom/ed  ( 34  kVV ) . 54 

Heine  2 1 _ C ' ompn rj s on  o I  !  I a  i ne  Spread . 56 

Figure  22  I’MMA  SuiTace  I cmperalure  (  nlculalcd  by  I l)S . 56 

I  inure  23  (  ompanson  oi  lolnl  Ileal  I  lu\  I  Vsinbiijjon  ;il  V  arious  Times . 57 

I'  leu  re  24 _ (.  '< * mmirison  of  M ; i ss  B u i  n  i ng  R  a t e . 60 

Heine  25  Comparison  Of  k . 64 

I inure  26  (.  omparison  ()|  a  . 64 


l; inure  27  Ileal  Release  Rale  (iivcn.Hy  I  I )S  (  aleulnlion  I  or  I  lie  Small-Scale  Test.. 65 
I  ieure  2S  Cross  (.  dmpai  ison  I  or  I  leal  Release  Rale  ( riven  By  I  DS  C  'aleulalion  l;or 

I  lie  l  aree-  Aiul  Small-Scale  I  csts . 66 

I  ieure  24  Cross  i  omparison  Between  Rredicied  I  Ian  to  Spread  Veloci  tv  By  SO  1  ;lli 


And  !•  I)S  And  I  xperimenlal  Vleasuremenls  By  t  )rloll  . 67 

I  ieure  5.0 _ I I  )S-(  alculaled  Surlaee  I  empcralure  Vs.  I  jeielil.  Varied  By  Time . 68 

He  tire  3  I  Incident  ( ‘onveeiive  I  leal  1  liry(  onijiarison,., . 69 


I' ieure  32 _ Incident  Rarlialive  I  leal  I  lu  \  (  pinjyirison  Between  I’rediclions  And  The 


l'Npenmcnliil  Values . 70 

I ’ieure  55  Numerical  Soluiion  ol  I  qualion  N5 . . . 75 


vi 


LIST  OF  TABLES 


l  ahlc  I  l-.xpenmcntal  Result*  and  I  lame  Height . 33 

I  able  2  I*  MM  A  Thermal  Properties  I  'sed  mjhe  I  l)S  Simulation . 44 

i  ah lo  3  Air  Thermal  Properties  I  ’sat  in  Ahmad  A  I  act h ' s  ( orrclalions . 59 

I  able  4  I  hermal  Properties  l  soil  In  S()l  IT  And  I  IDS . 63 


vii 
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1  INTRODUCTION 


Numerical  field  modeling  of  fires,  as  emerging  since  the  1980’s,  has  provided  an 
economic  and  reasonably  accurate  approach  to  the  assessment  of  fire  hazard.  At  the 
National  Institute  for  Standards  and  Technology  (NIST),  the  Fire  Dynamics  Simulator 
(FDS)  has  been  developed  based  on  the  fundamental  conservation  of  heat,  momentum 
and  mass  using  a  Large  Eddy  Simulation  (LES)  method.  The  effort  to  develop  FDS  and 
all  other  computer-based  models  has  been  paralleled  by  a  continual  effort  to  determine 
the  accuracy  of  the  predictions  by  comparing  the  models  to  laboratory  experiments.  As 
a  result  of  various  validation  studies,  the  FDS  has  proved  to  do  well,  especially  in  flow- 

i  i 

related  aspects,  i.e.  in  simulating  the  smoke  movement  and  smoke  plume  trajectory  . 
Spearpoint  et  al.3  simulated  a  fully  furnished  compartment  fire,  and  discovered 
unsatisfactory  predictions  in  his  flashover  scenario  of  a  room  fire. 


1.1  Motivation 

Fire  modeling  was  formally  introduced  to  the  fire  research  and  engineering 
arena  in  1960.  Along  with  the  introduction  of  fire  modeling,  the  knowledge  of  fire 
dynamics  has  grown  significantly.  This  knowledge  has  led  to  the  development  of 
complex  fire  models  that  can  only  be  practically  implemented  by  computers.  Many  of 
these  models  are  widely  used  in  Fire  Protection  Engineering  to  analyze  the  hazards  and 
risks  associated  with  building  fires  4.  Their  uses  have  become  widespread  in  the 
performance-based-design  applications. 


Fire  simulation  computer  codes  provide  a  means  to  evaluate  the  performance- 
based  design  for  fire  safety  under  a  range  of  possible  fire  scenarios.  With  wide 
applications  in  smoke  movement,  detector/sprinkler  activation,  fire  growth  and 
development,  to  building  structure  response,  the  fire  simulation  codes  are  helping  the 
fire  protection  engineering  arena  in  many  aspects. 

Jones  5  pointed  out  the  need  to  ensure  the  validity  of  fire  model  applications  in  a 
regulated  environment.  An  evaluation  process  for  model  validation  typically  includes 
the  following  components  6’7:  conceptual  validation,  model  verification,  and  operational 
evaluation.  Conceptual  validation  consists  of  determining  the  correctness  of  the 
assumptions  and  theoretical  basis  for  a  model  with  respect  to  its  intended  purpose. 
Model  verification  is  the  process  of  determining  the  correct  mathematical 
implementation  of  the  underlying  theory.  Operational  evaluation  is  quantifying  the 
accuracy  of  model  results.  The  first  two  processes,  conceptual  validation  and  model 
verification,  are  typically  performed  during  the  model’s  development  in  hope  of 
resolving  any  preliminary  issues  regarding  its  fundamental  validity.  The  last  process, 
operational  evaluation,  is  critical  to  establishing  the  acceptable  uses  and  limitations  of 
the  model. 

Evaluation  is  the  process  of  determining  the  models  range  of  application  and 
consists  of  comparing  model  output  to  reliable  data.  There  is  a  wide  range  of  possible 
scenarios  for  model  evaluation.  Hence,  this  process  has  historically  been  a  post¬ 
development  endeavor  that  continues  throughout  the  model’s  life. 
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In  a  recent  project,  Ma  31  explored  the  accuracy  and  limitations  of  FDS  on  axi- 
symmetric  fire  plumes.  They  examined  the  FDS  for  unconfined  fires  of  different 
dimensionless  heat  release  rate  Q*.  The  gridsize  plays  an  important  role  in  deciding  the 
flame  structure.  As  the  simulation  resolution  was  defined  as  the  grid  size  over  the 
characteristic  length,  they  found  that  the  optimum  resolution  of  a  pool  fire  simulation  is 
0.05.  In  other  words,  the  characteristic  length  is  defined  as 

2/5 

_ _ _  _  (n ’  |'5  .£) ,  and  the  grid  size  may  be  defined  as  5%  of  z*  so  as  to 

obtain  accurate  results.  Their  work  was  fundamental  to  further  work  on  the  flame- 
spread  prediction  by  FDS. 

o 

Lewis,  Rubini  and  Moss  performed  a  field  modeling  of  flame  spread.  Using 
the  SOFIE  code  9 ,  a  simplified  model  of  non-charring  solid  pyrolysis  has  been  coupled 
with  detailed  descriptions  of  turbulent  combustion,  soot  production  and  radiative 
exchange  to  create  a  field  model  for  fire  spread  prediction.  The  sensitivity  of  the 
predictions  to  alternative  soot  models  and  to  the  level  of  resolution  employed  is 
demonstrated  by  comparison  with  established  two-dimensional  experiments  on  PMMA 
slabs  reported  in  the  literature  by  Orloff  et  al  36.  The  SOFIE  Code  uses  a  two-equation 
k-'£  turbulence  model  for  modeling  of  buoyant  turbulent  flow,  while  the  FDS  uses  the 
LES  as  its  turbulence  model.  It  was  reported  that  the  level  of  agreement  between 
measurements  of  flame  spread  velocity,  mass  loss  rate  and  incident  heat  flux  is 
generally  good. 
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On  another  recent  related  research  by  Wang  et  al.  at  the  Universite  de  Poitiers  l0, 
where  non-spread  vertical  wall  fire  tests  with  a  narrow  range  of  Q*  are  examined,  it  was 
found  out  that  the  FDS-computed,  time-averaged  flame  height,  velocity  and  temperature 
profiles  are  relatively  in  good  agreement  with  the  measured  value  by  Orloff  ,  but  the 
deviation  increases  as  the  height  increases.  The  simulation  results  in  this  work  match 
their  findings. 

Spearpoint  et  al,  3  investigated  the  FDS  simulation  of  a  fully  furnished 
compartment  fire,  using  an  earlier  version  of  FDS,  and  compared  to  the  zone  model  Fire 
MD.  The  field  model  compared  unsatisfactory  with  zone  model  simulation,  partly 
because  of  inadequacy  of  the  sub-models  of  radiation  and  boundary  conditions. 

Floyd  11  performed  simulations  comparing  model  predictions  with  full-scale 
nuclear  facility  fire  test  data.  The  code  did  poorly  for  situations  where  the  fire  is  under¬ 
ventilated  or  in  the  flashover  stage.  The  weaknesses  were  attributed  to  the  crude 
radiation  and  combustion  models. 

A  series  of  simulations  were  performed  by  Friday  and  Mowrer  6  to  validate  FDS 
with  data  from  Factory  Mutual  and  Sandia  National  Laboratories  for  mechanically 
ventilated  enclosure  fires.  Functional  analysis  is  performed  for  comparing  results.  A 
grid  dependent  plume  study  was  performed  to  investigate  the  relationship  of  grid 
resolution  and  plume  temperature  predictions.  Increasing  the  number  of  grid  cells  in  the 
plume  region  does  not  improve  temperature  predictions  significantly,  but  does  result  in 
more  accurate  simulation  of  plume  turbulent  structures.  This  laid  a  solid  foundation  for 
our  grid-dependent  study  in  the  combusting  region. 
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In  the  newly  released  version  2.0  of  FD  S  debuted  on  December  4,200 1 ,  the  new 
features  mainly  consist  of  substantial  improvements  in  combustion,  radiation  transport 
and  the  boundary  conditions  . 

In  the  new  combustion  model,  a  more  comprehensive  method  that  handles 
oxygen  consumption  more  naturally  solves  an  equation  for  a  conserved  scalar  quantity, 
known  as  the  mixture  fraction,  which  is  defined  as  the  fraction  of  gas  at  a  given  point  in 
the  flow  field  that  originated  as  fuel.  The  model  assumes  that  combustion  is  mixing- 
controlled,  and  that  the  reaction  of  fuel  and  oxygen  is  infinitely  fast.  The  mass  fractions 
of  all  of  the  major  reactants  and  products  can  be  derived  from  the  mixture  fraction  by 
means  of  “state  relations,”  empirical  expressions  arrived  at  by  a  combination  of 
simplified  analysis  and  measurement. 

The  FDS  Version  2  handles  radiative  heat  transfer  in  a  better  way,  by  returning 
to  the  fundamental  radiation  transport  equation  for  a  non-scattering  gray  gas.  The 
equation  is  solved  using  a  technique  similar  to  finite  volume  methods  for  convective 
transport,  thus  the  name  given  to  it  is  the  Finite  Volume  Method  (FVM).  Using 
approximately  100  discrete  angles,  the  finite  volume  solver  requires  about  15  %  of  the 
total  CPU  time  of  a  calculation,  a  modest  cost  given  the  complexity  of  radiation  heat 
transfer. 

For  the  boundary  conditions,  the  few  changes  include  that  the  requirement  that 
all  internal  obstructions  (like  walls)  be  at  least  two  cells  thick  is  gone  in  FDS  2.  This  is 
due  to  a  revised  method  of  storing  information  about  solid  surfaces.  It  is  transparent  to 
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the  user,  but  does  result  in  less  “leakage”  of  heat  and  mass  through  one  cell  thick 


obstructions. 

However,  there  are  few  validations  and  the  subsequent  discussion  on  the  ignition 
and  flame  spread  on  FDS  predictions,  partly  because  the  earlier  versions  of  FDS  are 
more  smoke  movement  and  fire  plume  focused.  But  with  the  introduction  of  the  more 
advanced  sub-models  as  mentioned  above,  the  need  for  validation  and  further 
improvement  of  the  ignition  and  flame  spread  prediction  is  in  demand.  Much  of  the 
work  presented  here  focuses  on  the  ignition  and  flame  spread  issues,  plus  their 
relationship  with  the  incident  heat  flux  distribution. 


1 .2  Background 

For  a  computer  code  simulating  a  fire  scenario  with  the  emphasis  on  ignition 
and  flame  spread,  it  is  most  important  that  models  regarding  initiation  and  development 
of  the  burning  of  materials  are  properly  deployed.  Ignition  and  flame  spread  is  a  result 
of  the  thermal  decomposition  of  the  combustibles  and  the  ignition  of  the  resulting 
gaseous  pyrolysis  products  of  the  gaseous  mixture  above  the  combustible’s  surface.  It  is 
a  complex  phenomenon  that  involves  chemical  and  physical  processes  in  both  gas  and 
solid  phases.  The  combustible’s  thermal  properties  and  the  environmental  variables  are 
the  two  main  groups  of  factors  controlling  the  ignition  and  flame  spread  process.  The 
thermal  properties  usually  include  thermal  inertia  ( kpc ),  thermal  diffusivity  ( a  /.density 
(p)  or  surface  density  (pS),  heat  of  combustion  (ANC),  heat  of  vaporization  (A/-/v)  and 
pyrolysis  temperature  ( Tp)  etc.  The  environmental  properties  generally  include  incident 
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heat  fluxes,  geometry  factors,  flame  versus  gas-flow  orientation  etc.  While  the  incident 
heat  fluxes  consist  of  radiationhe-radiation,  convection  and  conduction  heat  flux,  thus 
further  making  the  scenario  a  more  complex  one.  Quintiere  provides  a  very  good 
description  of  this. 

Natural  fires  typically  involve  diffusion  flames  spreading  over  the  surface  of 
solid  combustibles.  The  spread  of  the  flames  results  from  the  complex  interaction  of 
transport  and  chemical  processes  that  occur  in  the  vicinity  of  the  boundary  separating 
the  burning  and  non-burning  regions.  For  a  flame  to  spread  over  the  surface  of  a 
combustible  material,  sufficient  heat  must  be  transferred  from  the  burning  region  of  the 
fuel  to  the  unhumt  material  to  cause  the  fuel  to  pyrolyze.  This  is  a  primary  condition 
that  must  be  met  for  the  flames  to  spread.  Once  the  fuel  pyrolyzes,  propagation  can 
procede.  The  gas  phase  reaction  of  fuel  vapor  and  oxidizer  introduces  two  central 
controlling  mechanisms  -  heat  transfer  to  the  unbumt  fuel  and  gas  phase  chemical 
reaction  -  are  often  separately  identifiable. 14 

In  the  fire  science  arena,  the  flame  spread  phenomena  is  broken  down  into  two 
stages  in  a  time  sequence:  ignition  and  flame  spread.  However,  the  controlling 
mechanism  mentioned  above  defines  whether  the  flame  is  self-sustainable  and  can 
propagate. 

The  gas  phase  ignition  from  a  solid  combustible  is  started  by  heating  up  the 
solid  by  an  incident  heat  flux,  which  can  he  a  heat  flux  from  the  flame  or  an  external 
source  independent  from  the  solid  combustible  in  question.  The  combustible  surface  is 
heated  up  and  raised  to  the  pyrolysis  temperature  ( Tp).  After  attaining  T„  the  pyrolyzed 
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vapor  leaves  the  surface,  is  diffused  and  convected  outwards,  mixes  with  the  oxygen, 
creating  a  flammable  mixture  near  the  solid  surface.  Flaming  ignition  will  occur,  if  the 
mixture  temperature  is  increased  so  that  the  energy  from  the  chemical  reaction  is  self- 
sustained. 

Once  the  fuel  is  ignited,  the  flame  propagates  across  the  surface,  establishing  a 
diffusion  flame  on  it.  The  flame  transfesr  heat  to  the  surface  and  more  fuel  is  vaporized, 
providing  the  necessary  gaseous  fuel  to  sustain  the  flame.  Three  independent  heat 
sources  lead  to  the  spread  of  a  flame:  gas  phase  heat  convection,  solid  conduction,  and 
external  heat  contribution.  Despite  the  heat  sources,  there  are  other  phenomena  which 
can  play  an  important  role  in  flame  spread,  such  as  air  flow  which  can  aid  or  delay  the 
flame  spread. 

Many  of  these  phenomena  have  been  considered  and  some  of  the  well- 
developed  theoretical  models  have  been  implemented  into  the  development  of  FDS, 
especially  Version  2. 


1.3  Approaches  and  Objectives 

In  this  work,  an  assessment  is  conducted  for  the  ability  of  the  FDS  model  to 
predict  heat  flux  to  a  vertical  wall  from  a  floor  burner  29 .  Furthermore,  FDS  prediction 
of  flame  spread  on  PMMA  is  examined  by  using  data  from  Factory  Mutual  33  as  the 
benchmark. 

A  static  case  of  no  spread  but  with  PMMA  burning  is  also  investigated.  It  is  to 
validate  the  dynamic  results  provided  by  FDS  with  the  steady  state  case  and  to  the 
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correlations  developed  by  Ahmad  and  Faeth  15  for  wall  fire  heat  flux  and  mass  burning 
flux. 

In  examining  flame  spread,  it  is  attempted  to  develop  a  simpler  model  for 
determining  the  surface  temperature  and  mass  flux  that  would  avoid  a  detailed 
computation  of  wall  conduction. 
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2  FIRE  DYNAMICS  SIMULATOR 


The  purpose  of  this  chapter  is  to  provide  a  description  of  the  fundamentals  for 
the  FDS  models  involving  ignition  and  flame  spread. 


2.1  General  and  Hydrodynamics  Model 

The  idea  of  numerically  simulating  a  fire  scenario  has  emerged  and  been 
continuously  developed  by  scientists  and  engineers.  Indeed,  the  fundamental 
conservation  equations  governing  fluid  dynamics,  heat  transfer,  and  combustion  were 
first  written  down  over  a  century  ago.  Despite  this,  practical  mathematical  models  of 
fire  (as  distinct  from  controlled  combustion)  are  relatively  recent  due  to  the  inherent 
complexity  of  the  problem.  The  governing  equations  are: 

Conservation  of  Mass: 

^P-  + V'pu  =  0  Equation  1 

at 

or:  ^tp(v-u)-0  Equation  2 

Conservation  of  Species: 

+  v  pT,u  =  V  ■  pDVY,  +  Wm  Equation  3 

at 

Conservation  of  Momentum: 

d(ftu)  + v  -puu  =  -  Vp  T pg  +  V •  T  Equation  4 

at 
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or:  p - =  -Vp  +  pg  +  V-r 


Equation  5 


Conservation  of  Energy: 

+  V'p/zu  =  —^-  +  q”'  +  V-kVT  Equation  6 

Bt  Dt 

Equation  of  State: 

p=DlpT  Equation  7 

When  rewriting  the  above  equations,  note  that  for  a  given  scalar  quantity  <p(x.y.z.t), 
because  of  the  mass  conservation  equation,  it  is  often  convenient  to  write  the  transport 
equations  between  these  two  forms: 

dp(p  Dd) 

— —  +  V  •  pu<p  =  p  — -  Equation  8 

dt  Dt 

While  dealing  with  the  governing  equations,  these  unknown  parameters  are 
solved  so  that  people  can  predict  a  fire  scenario.  These  to-be-solved  unknowns  are: 
density,/?;  velocity,  u;  enthalpy,  h  and  pressure,  p.  To  do  this,  the  following  parameters 
must  be  provided:  the  firepower  per  unit  volume,  q" ;  the  turbulent  viscous  stress,  T  ; 
and  the  thermal  conductivity  term,  V  iVf  . 


2.2  Diffusive  Terms  (LES) 

The  viscous  stress  tensor  in  the  momentum  equation  is  given  by: 


r 


=  p^def  u--j(v-u)l 


Equation  9 


11 


where  I  is  the  identity  matrix  and  the  deformation  tensor  is  defined  as: 


1  du_ 
dx 


1  /  du  dv  'j  1  /  du  dw  \ ' 

2  ^  dy  dx  J  2  ^  dz  dx  J 


def  u  ■  —  [vu  +  (Vu)  ]= 


]_ 

2 


' dv  dw  ^ 
Kdx  dy/ 


dv  1  /  dv  dw 

dy  2 1  dz  dy 


1  / dw  dw\  1  /  dw  dv\  dw 

2  ^  dJC  dz  J  2  d^  dz  j  dz 


Equation  10 


In  the  numerical  model,  there  are  two  options  for  treating  the  dynamic  viscosity 
p.  For  a  Large  Eddy  Simulation  (LES)  where  the  grid  resolution  is  not  tine  enough  to 
capture  the  mixing  processes  at  all  relevant  scales,  a  sub-grid  scale  model  for  the 
viscosity  is  applied.  Following  the  analysis  of  Smagorinsky  i6,  the  viscosity  can  be 
modeled  as: 


j 

-p(C,Aj^2(def  u)(def  u)-|(V-u)2 3jZ  Equation  11 

where  Cs  is  an  empirical  constant  value,  A  is  a  length  on  the  order  of  the  size  of  a  grid 


cell,  and  the  deformation  term  is  related  to  the  Dissipation  Function: 


s t  •  Vu  =  jitf  2(def  u)  (def  u)-^(V  u)r 


l 


du 


dx 


dv' 


2  —  +2  —  +2  — 


dy, 


dw\ 


dz 


dv  du) 

- h - 

dx  djy  / 


+ 


Equation  12 


+ 


/  dw  dw'' 
^ dz  dx  J 


2  ( du  dv  dw 

— - 1 - 1 - 

3  y  dx  dy  dz 


The  dissipation  function  is  the  rate  at  which  kinetic  energy  is  transferred  to 


thermal  energy.  It  is  a  source  term  in  the  energy  conservation  equation  that  is  usually 
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neglected  because  it  is  small  -  an  approximation  consistent  with  the  low  Mach  number 


equations. 

In  an  LES  calculation,  the  thermal  conductivity  and  material  diffusivity  are 
related  to  the  turbulent  viscosity  by 

^LES  =  /<::f  ;  (pd)/,les=~l  Equation  13 

Pr  Sc 

The  Prandtl  number,  Pr,  and  the  Schmidt  number,  Sc,  are  assumed  constants  for 
a  given  scenario.  There  have  been  numerous  refinements  of  the  original  Smagorinsky 
model  !7’1^19>  but  it  is  difficult  to  assess  the  improvements  offered  by  these  newer 
schemes.  There  are  two  reasons  for  this.  First,  the  structure  of  the  fire  plume  is  so 
dominated  by  the  large-scale  resolvable  eddies  that  even  a  constant  eddy  viscosity  gives 
results  almost  identical  to  those  obtained  using  the  Smagorinsky  model 20.  Second,  the 
lack  of  precision  in  most  large-scale  fire  test  data  makes  it  difficult  to  assess  the  relative 
accuracy  of  each  model.  The  Smagorinsky  model  with  the  empirical  fixed  value  Cs 
produces  satisfactory  results  for  most  large-scale  applications  where  boundary  layers 
are  not  well  resolved. 


2.3  Combustion  Models 

There  are  two  types  of  combustion  models  used  in  FDS.  The  choice  depends  on 
the  resolution  of  the  underlying  grid.  For  a  DNS  calculation  where  the  diffusion  of  fuel 
and  oxygen  can  be  modeled  directly,  a  global  one-step,  finite-rate  chemical  reaction  is 
most  appropriate.  However,  in  an  LES  calculation  where  the  grid  is  not  fine  enough  to 


13 


resolve  the  diffusion  of  fuel  and  oxygen,  a  mixture  fraction-based  combustion  model  is 
used.  In  this  work,  the  simulations  are  done  by  the  Mixture  Fraction  model. 

The  mixture  fraction  combustion  model  is  based  on  the  assumption  that  large- 
scale  convective  and  radiative  transport  phenomena  can  be  simulated  directly,  but 
physical  processes  occurring  at  small  length  and  time  scales  must  be  represented  in  an 
approximate  manner.  The  nature  of  the  approximations  employed  is  necessarily  a 
function  of  the  spatial  and  temporal  resolution  limits  of  the  computation,  as  well  as  the 
current,  often  limited,  understanding  of  the  phenomena  involved.  The  actual  chemical 
rate  processes  that  control  the  combustion  energy  release  are  often  unknown  in  tire 
scenarios.  Even  if  they  were  known,  the  spatial  and  temporal  resolution  limits  imposed 
by  both  present  and  foreseeable  computer  resources  places  a  detailed  description  of 
combustion  processes  beyond  reach.  Thus,  the  mixture  fraction  model  is  based  on  the 
assumption  that  the  combustion  is  mixing-controlled.  This  implies  that  all  species  of 
interest  can  be  described  in  terms  of  a  mixture  fraction  Z(\;  i).  The  mixture  fraction  is  a 
conserved  quantity  representing  the  fraction  of  material  at  a  given  point  that  originated 
as  fuel.  The  relations  between  the  mass  fraction  of  each  species  and  the  mixture  fraction 
are  known  as  “state  relations.”  The  state  relation  for  the  oxygen  mass  fraction  provides 
the  information  needed  to  calculate  the  local  oxygen  mass  consumption  rate.  The  form 
of  the  state  relation  that  emerges  from  classical  laminar  diffusion  flame  theory  is  a 
piecewise  linear  function.  This  leads  to  a  “flame  sheet”  model,  where  the  flame  is  a  two 
dimensional  surface  embedded  in  a  three  dimensional  space.  The  local  heat  release  rate 
is  computed  from  the  local  oxygen  consumption  rate  at  the  flame  surface,  assuming  that 
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the  heat  release  rate  is  directly  proportional  to  the  oxygen  consumption  rate,  which  is  an 
input  but  with  the  default  as  propane,  although  the  heat  release  rate  per  unit  mass  of 
oxygen,  AHo>  does  not  change  much  with  different  types  of  fuels. 

In  the  FDS  simulations,  the  combustion  heat  release  is  assumed  to  be  mixing 
controlled.  The  fuel  selected  for  model  validation  is  PMMA,  reflecting  the  wealth  of 
experimental  data  available  for  comparison  between  prediction  and  experiment.  PMMA 
is  assumed  to  thermally  pyrolyse  to  the  simple  monomer  CsHsCE  and,  at  its  simplest, 
chemical  reaction  then  proceeds  with  the  single  step: 

C5H)i02  +60,  =>  5C02  +  4H20  Equation  14 


The  mixture  fraction  Z  is  defined  as: 


Z  = 


sYf-(Yq-Y~) 

sy;  +  y~ 


Equation  1 5 


while  Z=1  in  the  fuel  stream  and  7j=Q  for  pure  air,  with  Y0  =  Y£  «  0.23,  and  for  PMMA 


monomer,  s  = 


_  (eX32)_1 
vfmf  OXioo)'  ' 
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Recalling  the  equations  of  specie  conservation  for  oxygen  and  fuel,  and  using 
the  following  transformations, 


sy;  +y; 


dyf 

Dt 


=  V  •  pDS'Yp  +  WF 


sy;  +  y; 


v .  PDv(r„  -  y;  V  w ; 


Equation  16 


Equation  17 


and  summing  up  the  above  two  equations,  yielding: 
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Equation  18 


D7 

p^-—^~  =  V  ■  pDVZ 
Dt 


The  assumption  that  the  chemistry  is  “fast”  means  that  the  reactions  that 
consume  fuel  and  oxidizer  occur  so  rapidly  that  the  fuel  and  oxidizer  cannot  co-exist. 
The  requirement  that  fuel  and  oxidizer  simultaneously  vanish  defines  a  flame  surface 
as:  (where  Z/-here  is  the  location  of  “flame  sheet”  or  stoichiometric  reaction) 


Z(xit)  =  Z,  and  Z  = 


r 


-  sy;  +  y* 

The  assumption  that  fuel  and  oxidizer  cannot  co-exist  leads 
relation”  between  the  oxygen  mass  fraction  Yo  and  mixture  fraction  Z, 


Vo(z)- 


vrrt-z/z,)  z<z, 


0 


z  >  z, 


Equation  19 

to  the  “state 


Equation  20 
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Figure  1  -Relationship  of  mixture  fraction  and  species  concentration  for 

burning  PMMA 
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The  oxygen  mass  conservation  equation  can  be  transformed  into  an  expression 


for  the  local  heat  release  rate  using  the  conservation  equation  for  the  mixture  fraction 
and  the  state  relation  for  oxygen  Yq(Z)\ 

P~~®'  -  d  dY0  DZ  _  pDVZ  ;  pDVY0  =  pD^VZ 


p  °  y  . 

Dt  dZ  Dt  dZ 


dZ 


Equation  2  1 


Yielding: 


-m-0.w(pD^WZ 


dYo 

dZ 


V  •  pDS/Z 


Equation  22 


The  main  purpose  of  mixture  fraction  model  is  to  obtain  the  heat  release  rate  per 
unit  volume,  q" ,  which  is  based  on  the  consumption  of  oxygen,  trig  : 


q"  _  A Ho  .  Equation  23 

The  mixture  fraction  model  has  provided  the  basis  for  the  combustion  model  by 
solving  only  one  additional  scalar  diffusion  equation,  and  the  infinite  rate  of  combustion 
is  consistent  with  LES  length  and  time  scales.  Furthermore,  all  the  product  species  can 
be  expressed  as  a  function  of  mixture  fraction. 

However,  the  assumptions  of  mixture  fraction  model  also  bring  some 
limitations.  Firstly,  the  diffusion  coefficient,  D,  has  to  be  same  for  all  gas  species. 
Secondly,  the  chemical  kinetics  is  simplified  to  one  step  reaction  with  infinite  reaction 
rate.  Thirdly,  fuel  and  oxygen  are  always  reacted  eventually,  regardless  of  temperature, 
thus  limiting  the  usage  of  this  model  to  certain  scenarios  -  such  as  the  ventilation- 
controlled  cases.  In  these  cases,  due  to  low  temperature  and  reduced  oxygen 
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concentration,  the  oxygen  and  fuel  does  not  react  as  predicted  by  the  Mixture  Fraction 
Model. 

2.4  Convective  Heat  Transfer  to  Walls 

The  heat  flux  to  a  solid  surface  consists  of  the  summation  of  gains  and  losses 
due  to  convection  and  radiation.  In  an  LES  calculation,  the  convective  heat  flux  to  the 
surface  is  obtained  from  a  combination  of  correlations  for  natural  and  forced 
convection: 

q"c  =  h-  A T  Equation  24 

q  v 

where  AT  is  the  difference  between  the  wall  and  the  gas  temperature,  C  is  the 
coefficient  for  natural  convection  (1.43  for  a  horizontal  surface  and  0.95  for  a  vertical 
surface)  ,  L  is  a  characteristic  distance  related  to  the  size  of  the  plate,  k  is  the  thermal 
conductivity  of  the  gas,  and  v  is  the  kinetic  viscosity  of  the  gas. 

2.5  Radiative  Heat  Transfer  to  Walls 

The  Radiative  Transport  Equation  (RTE)  for  a  non-scattering  gray  gas  is: 

s-V/A(x,s)=fc(x,Aj/,(x,s)-/(xJs)]  Equation  26 


i 

h  =  max  cIatIj ,0.037 
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where  I ,  (x,s)  is  the  radiant  intensity  at  wave  length  A,  lb  (x,s)  is  the  source  term  given 

by  the  Planck  function,  s  is  the  unit  normal  direction  vector  and  K(x)  is  the  gray-gas 
absorption  coefficient.  In  practical  simulations,  the  spectral  dependence  cannot  be 
solved  accurately.  Instead,  the  radiation  spectrum  is  divided  into  a  relatively  small 
number,  n,  of  bands,  and  a  separate  RTE  is  derived  for  each  band. 

s-V/„(x>s)=  r:„(x][4„(x,s)-/(x,s)J  Equation  27 

The  limits  of  the  bands  are  selected  to  give  an  accurate  representation  of  the 
most  important  radiation  bands  of  CO2  and  water.  When  the  intensities  corresponding  to 
the  bands  are  known,  the  total  intensity  is  calculated  by  summing  over  all  the  bands: 

N 

/(x,s)=  ^/„(x,s)  Equation  28 

Even  with  a  reasonably  small  number  of  bands,  the  solution  of  ARTE’s  is  very 
time  consuming.  Fortunately,  in  most  large-scale  fire  scenarios  soot  is  the  most 
important  combustion  product  controlling  the  thermal  radiation  from  the  fire  and  hot 
smoke.  As  the  radiation  spectrum  of  soot  is  continuous,  it  is  possible  to  assume  that  the 
gas  behaves  as  a  gray  medium.  The  spectral  dependence  is  lumped  into  one  absorption 
coefficient  (  N  =  1 )  and  the  source  term  is  given  by  the  blackbody  radiation  intensity: 

I h  (x  )  =  Equation  29 

JZ 

In  optically  thin  flames,  where  the  amount  of  soot  is  small  compared  to  the 

amount  of  CO2  and  water,  the  gray  gas  assumption  may  produce  significant  over- 

1* 

predictions  of  the  emitted  radiation  .  For  the  calculation  of  the  gray  or  band-mean 
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absorption  coefficients,  k>,5  a  narrow-band  model,  RADCAL/3,  has  been  implemented  in 
FDS.  At  the  start  of  a  simulation,  the  absorption  coefficients  are  tabulated  as  a  function 
of  mixture  fraction  and  temperature.  During  the  simulation,  the  local  absorption 
coefficient  is  found  by  table-lookup.  In  calculations  of  limited  spatial  resolution,  the 
source  term,  lh ,  in  the  RTE  requires  special  treatment  in  the  neighborhood  of  the  flame 

sheet  because  the  temperature  values  are  smeared  out  over  a  grid  cell  and  are  thus 
considerably  lower  than  one  would  expect  in  a  diffusion  flame.  Because  of  its 
dependence  on  the  temperature  raised  to  the  fourth  power,  the  source  term  must  be 
modeled  in  those  grid  cells  cut  by  the  flame  sheet.  Elsewhere,  there  is  greater 
confidence  in  the  computed  temperature,  and  the  source  term  can  assume  its  ideal  value 
there: 


koT 4 

Outside  flame  zone 

JT 

.  tn 

Equation  30 

XrQ 

Aji 

Inside  flame  zone 

Here,  q”‘  is  the  chemical  heat  release  rate  per  unit  volume  and  %r  is  the  local 
fraction  of  that  energy  emitted  as  thermal  radiation.  There  is  a  difference  between  the 
prescriptions  of  a  local  Xr  and  its  resulting  global  equivalent.  For  a  fire  with  D<  1  m,  the 
local  X/'is  approximately  equal  to  its  global  counterpart.  However,  as  the  fire  increases 
in  size,  the  global  value  will  typically  decrease  due  to  a  net  re-absorption  of  the  thermal 
radiation  by  the  increasing  smoke  mantle. 
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The  boundary  condition  for  the  radiation  intensity  leaving  a  gray  diffuse  wall  is 
given  as  the  following,  and  the  radiative  heat  flux  at  the  surface  is  obtained  from  the 
boundary  condition, 

1  -  e 

TV 

where  /H,(s)  is  the  intensity  at  the  wall,  e  is  the  emissivity  of  the  wall,  and  hw  is  the  black 
body  intensity  at  the  wall. 

The  radiant  heat  flux  vector  qr  is  defined  as: 

qr(x)  =  Equation  32 

The  radiative  loss  term  in  the  energy  equation  is 

-V  .q,(x)=fc(x)j£  /(x,s)/Q -4a7,(x)t  Equation  33 

In  words,  the  net  radiant  energy  gained  by  a  grid  cell  is  the  difference  between 
that  which  is  absorbed  and  that  which  is  emitted. 


Equation  3  1 


4(s)=^  + 


2.6  Pyrolysis  Model  -  Thermally  Thick  Solid 

This  work  only  discusses  the  pyrolysis  model  for  a  thermally  thick  solid,  which 
is  also  referred  to  semi-infinite  solid.  For  a  thermally  thick  solid,  the  slab  is  heated  on 
one  side  only,  with  insignificant  heat  losses  from  the  rear  face.  In  many  fire  engineering 
problems  involving  transient  surface  heating,  it  is  adequate  to  assume  "thermally  thick 
behavior"  if  A>6,  where  A  is  the  thickness  of  the  solid  in  interest,  6  is  the  thermal 
penetration  depth  and  5  a  -Jctf  24 . 
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Firstly,  the  solid  is  heated  by  incident  heat  flux  and  the  surface  temperature  is 
increased.  One-dimensional  heat  conduction  equation  for  the  material  temperature,  TJx; 
t)  is  applied  in  the  directionx  pointing  into  the  air/solid  interface  (x  -  0): 

p  cK  —  =  ks  ~~~r  -  ■  •  for  T  <  Tp  Equation  34 

dt  dx* 

c)  T 

-  ks  =  q "  +  q”  -  q"RR  Equation  35 

OX  ,mrfuce 

where  ps ,  cs  and  ks  are  the  constant  density,  specific  heat  and  conductivity  of  the  solid; 
q" ,  q”  and  q"RR  are  the  convective,  radiative  and  re-radiation  heat  fluxes  at  the  surface. 

While  the  heat-up  proceeds  to  a  level  that  the  surface  temperature  attains  to  a 
prescribed  pyrolysis  temperature  of  the  solid,  Tin  the  boundary  conditions  of  the  above  2 
equations  change  to: 

3  J1 

Ts (0,0=7;,  q"p  =  q"  +  q’r  -  qRR  ~  ks  Equation  36 

OX  surface 

Here,  q"p  is  energy  available  for  pyrolyzing  fuel  via  the  simple  relationship: 

•  n 

m"  =  ^  —  Equation  37 

A Hv  4 

where  m"  is  the  burning  rate  or  mass  loss  rate,  while  A//,,,  is  the  heat  of  vaporization, 
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2.7  Proposed  Algorithm  for  Surface  Temperature  of  Thermally-thick  Material 


In  FDS,  a  solid  has  thickness  (5 .  This  thickness  is  divided  into  N  cells. 
Temperature  in  ith  cell,  Tsj,  is  obtained  in  time  with  a  Crank-Nicholson  scheme: 


rprl  +  l  rrtfl 

1  s,i  J  sj 

6t 


a 


( T"i+l  -  277  +  77, 


.v,/  1 


+t,%\ 


& 


Equation  38 


where  Is  /  s  N  .  Boundary  condition  is  discretized  as: 

T  -T 

-  — —  =  q'l  +  q"  -  q"RR  Equation  39 

to* 

and  wall  temperature  Tw  =  (  TSj0  +  Ts,i  )  /  2. 

As  5  is  defined  by  the  user,  and  N  is  20  by  default,  the  accuracy  of  the  results  by 
the  Crank-Nicholson  scheme  depend  upon  the  resolution  of  the  FDS  grid  size. 

An  alternative  algorithm  is  proposed,  that  has  no  dependency  upon  the 
resolution  of  grid  size.  It  is  based  on  the  same  physical  model  of  1-D  heat  conduction, 
i.e.  equation  29.  A  solution  was  developed  to  estimate  the  surface  temperature  when  a 
solid  is  exposed  to  a  time-related  net  surface  heat  flux  q"er(v)=  q"c  +  q"r  -q"RR  : 


Ts  -  Tx  ~  ,  Equation  40 

-qjikpc^^t-r 

In  order  to  solve  the  above  equation,  the  incident  heat  flux  as  a  function  of  time 
has  to  be  provided  as  a  given  input. 

The  detailed  derivation  of  the  solution  for  the  above  equation  is  enclosed  as 
Appendix  A,  and  the  solution  is: 


23 


Comparing  to  the  existing  algorithm  for  surface  temperature,  the  proposed 
algorithm  is  independent  of  the  resolution  of  FDS  domain.  Furthermore,  this  is 
computationally  cost-effective. 


2.8  Proposed  Algorithm  for  Burning  Rate  of  Thermally  Thick  Material 

25  26  27 

The  following  algorithm  was  developed  and  presented  by  Quintiere  et  al.  ’  ' 
using  a  one-dimensional  heat  transfer  model. 

First,  the  preheating  period  before  pyrolysis  is  examined  to  obtain  the  thermal 
penetration  depth  which  is  a  function  of  time.  Again,  the  one  dimensional  unsteady 
conduction  equation,  as  Equation  34,  applies;  plus  the  constant  initial  temperature, 

t  =  0,  T  =  Tx  Equation  42 

and  at  the  surface. 

-  k  s  q”«,  =  [-  K  fa  -  Tg  )}-  q"f,  +  (l  -  r fc  -  eoTs 4  Equation  43 

,  / solid, y-0 
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Employing  a  standard  approximate  integral  solution  to  solve  the  above  problem, 
at  some  penetration  depth,  6,  from  the  surface,  the  temperature  is  and  the  heat  flux  is 


zero,  i.e.: 


or 


y  =  5,T  =  Toand  —  =0 
fly 


Equation  44 

The  integral  equation  is  derived  by  integrating  Equation  34  for  0^  y  £  <5  (t ) .  It 


can  be  shown  that  it  becomes: 


dt 


■jffr-r.)*--* 


ldT\ 


Equation  45 


y-0 


A  polynomial  profile  for  T  within  8  is  assumed: 


T  =  Ay2  +  By  +  C 


Equation  46 


A,  B,  and  C  are  determined  by  applying  the  boundary  conditions  given  by 
Equation  42,  Equation  43,  and  Equation  44.  This  follows  as: 


A  =  -~ 
26 


C  =  T  - 


B_ 
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kB  +  hc^d+  q"f'r  +  ^  "  r^"  '  ea(r-  “  =  0 


Equation  47 


Equation  48 


Equation  49 


Substituting  this  profile  into  Equation  45,  along  with  the  initial  conditions: 


t=0,  8=0 


Equation  50 


yields, 
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5  =  V6 at 


Equation  51 


which  is  independent  of  the  parameter  B. 

The  complete  solution  is  found  by  solving  Equation  49  at  each  time  by  a 
Newton-Raphson  method  to  find  B.  Then  the  solution  is  carried  out  in  time  until 
T(0,t)=C=Tp.  The  values  of  6  and  the  temperature  profile  at  the  time  of  the  onset  of 
pyrolysis  provide  the  initial  conditions  for  the  gasification  period  discussed  next. 
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Figure  2  -  Heat  &  Mass  Transfer  Processes  for  a  Thermally  Thick  Material 

Then,  the  gasification  period  following  the  preheating  stage  is  discussed.  The 
governing  equations  for  the  gasification  period  of  thermally  thick  material  can  be 
derived  in  integral  form  from  the  control  volumes  (CVI  and  CV2)  selected  in  Figure  2 
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for  the  vaporization  plane  and  solid  respectively.  The  integral  equation  can  also  be 
derived  by  integrating  the  partial  differential  equation  of  Steckler  et  al.  28  over  the 
penetration  depth,  6. 

Conservation  of  mass  applied  to  the  solid  CV2  expresses  that  the  rate  of  mass 
changed  in  CV2,  which  is  defined  by  an  increasing  6,  is  equal  to  the  rate  of  mass 
entering  the  control  volume  at  y  =  6.  Conservation  of  energy  applied  to  CV2  considers 
enthalpy  transport  across  the  moving  vaporization  surface,  before  vaporization  occurs, 
by  the  mass  flux  pv  (v  is  the  surface  regression  velocity,  v=dys/dt);  and  enthalpy 
transfer  across  the  surface  defined  by  5  due  to  the  mass  flux,  p  (dS/dt  +v  ),  which  enters 
the  CV.  Specific  heat  for  the  solid  at  constant  pressure  and  volume  are  assumed  equal 
and  the  reference  state  for  enthalpy  is  taken  at  Tx  .  Hence,  for  CV2: 


solid,y-0 


Conservation  of  energy  for  CVI  yields: 


With  the  boundary  conditions: 

y  =  0,  T  =  Tv 

y  =  6,  T  -  T0  and  —  =0 

dy 


Equation  52 


h "AH y  =  [- h,  —  j  +  q)r  +  (l  - r)q"r  - soT*  Equation  53 

\  ^  /  solid, y-0 


Equation  54 


Again,  a  polynomial  temperature  profile  within  6  is  assumed  for  the  integral 
solution: 
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T  =  Ay2  +  By  +  C 


Equation  55 


Using  the  three  boundary  conditions  in  Equation  54,  the  parameters  A,  B,  and  C 
are  obtained: 


4  =  £=-2^-rJ,  C-T 

o  6 


Equation  56 


Substituting  Equation  55  and  Equation  56  into  Equation  52,  it  becomes: 

pc^r-T.)dd  ,,Cr  2 

- L. - -  +  m  c\f  -T„)= - : - .or. 


dt 


1  dz.  +^L_  =  2a_ 
3dt  p  a 


Equation  57 


The  initial  conditions  at  the  onset  of  pyrolysis  are  the  values  when  the  surface 
temperature  of  the  heat-up  process  achieves  T,.  And  the  initial  value  of  6  in  Figure  2 


right  before  pyrolysis  is  6((  =  t,  P  .  In  summary,  the  initial  conditions  at  the 


onset  of  pyrolysis  are: 

t=t„  Ts=Tpand<$(=0=V^ 

Let q"net  (t  ) s  q"  +  q”  - q”RR  as  the  net  surface  heat  flux,  or: 

r„,  4 -A, (7;  -T  +(1 


Equation  58 


Equation  59 


Substituting  Equation  59  into  Equation  53,  and  substituting  Equation  56  into  the 
3  T 

derivative  of  —  ,  Equation  53  can  be  rewritten  as: 

dy 


28 


/  ,dT\  .„  2k{fp-Tx) 

=qnel-m  Aff„  = - - - 

l  ay  ,a/M  d 


Equation  60 


Explicitly, 


*?»<■/ 


2k(Tp-T.) 


m  - 


Equation  61 


A Hv 


Utilizing  the  initial  condition  at  pyrolysis  as  Equation  58,  and  discretizing 
Equation  57,  the  burning  rate  algorithm  is  summarized  as: 


Before  ignition:  t  £  tt  : 


+  !  _  rp 

*  p  ■*  QC  * 


1 


t]jt  ■  kpc 


m"  =  0 


(a  +  h ) 


Equation  62 


After  ignition:  t  >  tjg 


^n  +  l  _ 


T.  =T, 

( 6 a  3m 


i 

■  a 


1  = 


A Hr 


P 

nei,n  +  \ 


■M  +  S, 

2*(r„  -r.  > 


Equation  63 


In  summary,  the  existing  FDS  algorithm  utilizes  similar  physics  while  the 
proposed  model  eliminates  the  need  for  a  repeated  solution  of  the  PDE.  A  benefit  is 
reduced  computation  cost.  The  potential  weakness  of  the  proposed  model  for  mass 
burning  rate  is  that  it  is  for  a  non-charring  semi-infinite  solid. 

The  proposed  algorithm  is  not  incorporated  into  the  FDS  code  and  is  a  subject  of 
future  study  for  quantifying  its  computational  benefits,  and  the  following  evaluations  on 
the  FDS  code  are  based  on  the  existing  models. 
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3  BACK’S  WALL  FIRE  TESTS 


2o 

The  series  of  experiments  used  in  the  study  was  performed  by  Back  et  al. 
These  fire  tests  showed  that  the  peak  wall  heat  fluxes,  for  fires  of  50-500  kW,  could  be 
40-120  kW/m2.  Such  heat  fluxes  are  sufficient  to  ignite  most  combustible  materials.  It 
is  crucial  that  the  heat  flux  prediction  by  FDS  can  be  maintained  at  a  reasonably 
accurate  level,  in  order  to  achieve  reasonable  results  for  surface  flame  spread  and 
subsequent  flashover  prediction.  This  work  is  to  provide  a  comparison  between  the  FDS 
prediction  and  the  experimental  data. 


3.1  Description  of  the  Wall  Fire  Experiment 

Incident  wall  heat  flux  distributions  have  been  measured  for  square  propane 
burner  fire  sources  placed  against  the  wall.  Peak  heat  fluxes  of  40-120  kW/m2  were 

'y 

measured  for  energy  release  rates  of  50-500  kW/m  and  burner  edge  lengths  of 
0.28-0.7  m.  Peak  heat  flux  was  found  to  be  a  strong  function  of  energy  release  rate  hut 
were  insensitive  to  burner  size.  The  experimental  apparatus  is  shown  in  Figure  3.  The 
specific  test  conditions  are  included  in  Table  1. 
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Thermocouple  1.25 
cm  from  wall  surface, 
heat  flux  meter  on 
surface;  15  cm  apart. 


1.25  cm  t  h  i  c  k 
gypsum  board  back 
by  steel  frame 


Square  propane 
burner  w/  0=50-500 
kW,  edge  0.28-0.7  m 


Propane  supply 


Figure  3  -  Experimental  Setup  for  Beck's  Wall  Fire 


Eight  Medtherm  total  heat  flux  transducers  were  installed  flush  with  the  wall 
surface.  A  vertical  array  of  thermocouples  was  installed  beside  the  total  heat  flux  array, 
approximately  1 .25  cm  from  the  wall  surface. 


3.2  Back's  Experimental  Results 
Flame  Height 

Flame  heights  were  measured  by  visual  analysis  of  the  videotape  and 
thermocouple  measurements.  The  time-averaged  and  50%  visual  intermittence  results  of 
these  measurements  are  shown  in  Table  1 . 
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Centerline  Temperatures 

After  adjusting  the  data  for  the  experimental  virtual  source  (zo),  the  data  appear 
to  lie  in  two  discrete  regions,  the  flame  and  the  plume.  The  temperatures  in  the 
continuous  combustion  region  are  approximately  constant  (900  “C).  The  temperatures  in 
the  non-reactive  zone  decrease  as  a  function  of  height  and  the  energy  release  and  are 
well  represented  by  the  following  correlation: 

ATcl  =40.6Q2i(z-z0)5n  ...for  (z  -  z(t)j  Q"  5  >  0,16  Equation  64 

Beck  et  al.  also  found  out  that  their  experimental  virtual  source,  zo,  was  in 
agreement  with  Heskestad’s  correlations: 

z0  =  0.083Q"75  -  1 .020  Equation  65 

Heat  Flux 

The  centerline  measured  heat  flux  profiles  are  shown  in  Figure  2  with  the  height 
normalized  by  the  flame  height.  The  flux  is  approximately  equal  to  the  peak  flux  over 
the  lower  40%  of  the  flame  height  and  decays  to  20  kW/m2  at  the  flame  height.  The 
consistency  of  the  20  kW/m2  incident  flux  at  the  flame  height  and  its  correspondence 
with  a  blackbody  at  500  °C  is  good,  as  500  °C  matches  the  visible  flame. 
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ZJLf 

Figure  4  -  Measured  centerline  heat  flux 


Test 

Number 

Q 

(kW) 

Z* 

(m) 

& 

50%  Visual 
Intermittence 
Flame  Height  (m) 

Time-averaged  | 
500  °C 

Flame  Height  (m) 

1 

0.298 

0.867 

0.79 

0.78 

2 

56 

0.70 

0.305 

0.093 

0.36 

0.33 

3 

68 

0.48 

0.330 

0.51 

4 

106 

0.37 

0.394 

0.864 

1.00 

1.02 

5 

136 

0.48 

0.435 

0.578 

0.87 

0.86 

6 

204 

0.48 

0.512 

0.868 

1.45 

1.45 

7 

220 

0.70 

0.527 

0.364 

1.20 

1.29 

8 

313 

0.866 

2.20 

2.29 

9 

523 

0.746 

0.866 

— 

2.90 

Table  1  -  Experimental  Results  and  Flame  Height 
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4  EVALUATION  OF  THE  FDS-PREDICTION  ON  BACK  S  WALL 


TEST 


4.1  FDS  Input 

The  geometry  data,  i.e.  the  wall  and  the  burner,  was  taken  directly  from  the 
experimental  setup.  The  domain  consists  of  a  solid  wall  against  the  burner  on  the  floor, 
and  all  other  walls  are  set  as  openings  to  the  ambient.  Ma  has  recommended  that  the  cell 
size  be  5%  of  z*,  the  characteristic  plume  length  scale,  for  accurate  simulation  results 
using  the  FDS  31 .  And  z*  is  represented  by: 


Q 


ter- 


D 


Where  in  this  case,  the  ambient  temperature  is  selected  as  20  °C. 


Equation  66 


Figure  5  -  FDS  Simulation  Domain  for  Back's  Wall  Fire  Test 
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In  the  test,  the  wall  was  constructed  of  gypsum  wallboard  with  a  very  low 


conductivity  and  back-surface  heat  loss.  Keeping  this  in  mind,  a  wall  with  actual 
gypsum  thermal  properties  was  used  with  the  values30  of  k=0.48  W/m-K,  _=4. 1X10'7 
m2/s  andthickness=0.013  m.  This  unreleased  version  of  FDS  used  in  this  paper,  V.  1.99 
dated  April  30th  2001,  has  incorporated  variable  specific  heat  as  a  function  of 
temperature. 


42  FDS  Predictions  on  the  Centerline  Temperature 

The  time-averaged  centerline  temperature  was  correlated  by  Equation  64.  That 
was  measured  1.25  cm  from  the  wall,  at  the  center  of  the  burner.  The  corresponding 
results  were  computed  from  the  FDS.  Runs  for  each  test  configuration  in  Table  1  were 
computed  for  10  seconds  of  simulation  time.  This  was  done  on  a  Pentium3-866MHz, 
384MB  computer  which  took  1-1.5  hour  for  about  each  simulation  second.  The 
corresponding  time-averaged  temperature  were  computed  from  3  to  10  seconds  of 
simulation  time.  The  comparison  between  the  experimental  correlation  and  the  FDS 
simulation  is  shown  in  Figure  6,  the  numbered  curves  correspond  to  the  Test  Number  in 
Table  1 .  Both  sets  descend  while  Q"d  decreases,  hi  the  lower  part  of  the  flame  region, 
the  model  over-predicts  by  20-30%  the  gas  temperature,  which  is  empirically  defined 
by  a  fixed  temperature  of  900  °C,  hi  the  upper  part  of  the  flame  region  and  in  the  plume 
region,  the  FDS-predicted  temperature  falls  quickly  and  the  result  is  poor. 
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4.3  FDS  Predictions  on  the  Wall  Heat  Flux  Distribution 

As  shown  in  Figure  7,  the  model  successfully  predicts  the  maximum  flux,  which 
is  located  just  above  the  fuel  surface.  However,  similar  to  the  predictions  for 
temperature,  the  predictions  for  heat  flux  level  decay  quickly.  For  example,  the 
experimental  data  showed  that  flame  height  is  at  1.45  m  for  Test  6;  however,  at  this 
height,  the  model  predicts  a  total  heat  flux  as  low  as  5  kW/m__,  which  is  substantially 
lower  than  the  measured  total  heat  flux  of  20  kW/m_  at  the  flame  height. 

As  shown  in  Figure  8,  the  heat  flux  was  broken  down  as  radiative  heat  flux  and 
convective  heat  flux,  and  compared  to  the  experimental  results. 
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The  distribution  of  radiation  heat  flux  and  convective  heat  flux  is  satisfactory. 
Convective  heat  fluxes  are  predicted  at  maximum  20  kW/rrC  at  the  base  of  the  wall,  and 
it  gradually  decreases  to  2-5  kW/m2  while  the  height  increases.  The  radiative  portion, 
that  is  dominant,  is  decreasing  quickly  as  the  height  increases,  and  this  finding  may 
suggest  further  improvement  in  the  radiation  sub-model  he  needed. 


Figure  7  -  Heat  Flux  Distribution 
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Z  (m) 

Figure  8  -  Heat  Flux  Distribution  for  Test  6 


4.4  FDS  Predictions  on  the  Flame  Height 

Figure  7  plots  the  dimensionless  results  for  flame  height.  The  data  are  indicated 
by  the  test  numbers.  The  correlation  for  the  axisymmetric  freestanding  fires  is  given  by 
Heskestad: 

zf  =  0.23Q215  - 1 .02D  Equation  67 

The  simulation  flame  height  was  based  on  the  time  averaged  from  3  to  10 
seconds  for  the  locus  of  the  energy  released  rate  at  which  it  achieves  99,  99.9,  and 
99.99%  of  its  full  value.  Ma  31  found  that  the  99.99%  locus  best  matched  axisymmetric 
free  tire  plume  data. 
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4.5  Summary  of  the  FDS  Predictions  on  Beck's  Wall  Fire  Tests 

The  FDS/LES  model  gives  results  for  temperature,  wall  heat  flux,  and  flame  in 
qualitative  agreement  with  data.  The  temperature  and  the  heat  flux  are  within  about 
30%  of  measured  results  for  peak  values  near  the  base  of  the  wall,  but  both  drop  much 
quicker  than  the  experimental  results.  The  weakness  appears  to  be  combustion  model 
which  is  giving  too  high  temperature  predictions  inside  the  flame  region  (Figure  6). 
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5  FM’S  LARGE  SCALE  PMMA  WALL  FIRE  TEST 


The  benchmark  data  set  to  be  compared  with  the  FDS  prediction  was  developed 

in 

by  Factory  Mutual  Research  Corporation.  Tewarson  performed  a  comprehensive 
study  on  the  flammability,  flame  spread,  and  flame  extinction  behaviors  of 
polymethylmethacrylate  (PMMA)  slabs.  A  couple  of  years  later,  Wu  et  al.  33  at  FMRC 
continued  the  research.  In  validating  FMRC’s  Flame  Spread  &  Growth  model,  as  well 
as  providing  badly  needed  full-scale  heat  flux  data,  an  upward  fire  spread  experiment 
was  carried  out  under  the  Fire  Products  Collector  at  the  FMRC  Test  Center. 


5.1  FMRC  500-kW-Scale  Flammability  Apparatus  and  the  Experimental  Setup 

In  the  experiments  32,  polycast  PMMA  with  a  density  of  1190  kg/m3  was  used. 
The  500-kW  Scale  Flammability  Apparatus  consists  of  a  lower  section  and  an  upper 
section. 

The  lower  section  is  used  for  the  measurements  of  time  to  ignition,  mass  loss 
rate  during  pyrolysis  and  combustion,  pyrolysis  and  flame  heights,  and  flame  spread 
rate,  as  well  as  concentrations  of  Yo  and  gaseous  agents  and  water  application  rate 
required  for  flame  extinction.  The  upper  section  is  used  for  the  measurements  of  total 
mass  and  volumetric  flow  rates  of  chemical  compounds  -  air  mixtures,  temperature, 
concentrations  of  CO,  CO2,  02,  total  gaseous  hydrocarbons,  HjO,  particulates,  and  other 
compounds  such  as  HCI,  NO,  and  N02,  and  optical  transmission,  corrosive  nature  of  the 
chemical  compounds-air  mixture,  and  electric  charge  on  the  particulates. 
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The  full-scale,  upward  fire  spread  experiment  was  conducted,  using  a  0.025  m 
thick  PMMA  wall,  0.58  m  wide  x  5  m  high.  The  PMMA  wall  was  extended  another  0.3 
m  on  each  side  by  Marinite  panels.  At  the  outer  edge  of  the  Marinite  panels,  a 
perpendicular  0.6  m  flow  barrier  (24  gauge)  steel  is  used  to  minimize  the  effects  of 
room  drafts.  On  the  top  of  the  wall,  a  3  m  extension  provides  a  background  for 
measuring  the  flame  heights.  Among  many  data  acquisition  instruments,  seven  total 
heat  flux  gauges  and  thermocouples  were  placed  at  various  on  the  PMMA  wall. 


To  Fire  Product  Collector 


1.2  m 


Flame  Extension 
Board  (24  gauge 


steel) 


PMMA 


Marinite 


Flow  holder 

(24  gauge  steeL  ^ 

0.6  m  wide) 


Load  Cell 


3m 


5  m 


Figure  1 0  -  Full  Scale  PMMA  Wall  Fire  Test 33 
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Symbol 

Equipment  List 

1 

Blast  gate  and  linear  actuator 
control 

2 

Particular  sampling;  pressure  and 
corrosion  measurement 

3 

Optica1  transmission  measurement 

4 

Gas  temperature  measurement 

5 

Vertical  testing  sample  slab 

6 

Infrared  heaters 

7 

Load  cell 

8 

Air  and  oxygen  inlet 

Figure  1 1  —  FMRC  50-kW  Flammability  Apparatus  as  shown.  The  500-kW  apparatus 

similar  except  for  size. 32 


The  experiment  was  performed  in  the  500-kW  Scale  Flammability  Apparatus, 

using  vertical  slabs.  To  hold  the  samples  in  a  vertical  orientation,  a  ladder-like  sample 

holder  was  used,  as  shown  in  Figure  1 1 .  Flame  spread  was  restricted  to  a  single  surface 

by  tightly  covering  the  back  of  the  slab  with  a  3-mm  thick  ceramic  paper. 

For  ignition  of  flame  spread  with  a  uniform  pyrolysis  front  for  each  test,  a 

0.025-m-high,  0.025-m-wide,  0.10-m-long  aluminum  dish  containing  about  10  ml 

methanol  was  placed  at  the  bottom  of  the  slab  and  ignited  by  a  match. 

For  the  measurement  of  pyrolysis  height  and  flame  spread  rate,  lines  were 

marked  at  intervals  of  25  mm  on  the  surface  across  the  width  of  the  slabs.  For  the 

measurement  of  flame  height,  lines  were  marked  every  25  mm  on  the  sample  holder  and 

holes  were  made  every  25  mm  on  the  aluminum  extension.  A  stopwatch  was  used  to 

record  visually  times  for  pyrolysis  front  and  flame  tip  arrival  at  each  line. 

Most  of  the  flame  spread  experiments  were  performed  at  vg=  0.09  m/s  with  vg  as 

the  environment  air  velocity;  the  other  value  used  was  0.18  m/s.  These  values  were  not 

expected  to  affect  the  flame  spread  rate,  as  they  are  significantly  less  than  ly  =  0.30  m/s, 

1 1 

at  or  above  which  flame  spread  is  affected  by  the  gas  flow  . 

In  the  experiments,  mass  loss  rate  during  flame  spread  was  not  measured. 
However,  initial  and  final  weights  of  the  samples  were  recorded. 

In  the  upper  section,  all  the  compounds  generated  during  pyrolysis,  combustion, 
flame  spread  and  extinction  were  captured,  along  with  ambient  air  in  the  sampling  duct. 
The  upper  section  operates  as  large-scale  cone  calorimeter,  so  this  is  not  of  the  interest 
of  this  work. 
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Thermal  Properties 

Value  and  Unit 

k,  conductivity 

0.209  W/m-K 

cp,  specific  heat  of  the  solid 

1.46kJ/kg-K 

p,  density  of  the  solid 

1 190 kg/m 3 

kpCp,  thermal  inertia 

0.364  kW2-s/m4-k2 

a,  thermal  dijfiusivity 

1.2E-7m_/s 

Tp,  pyrolysis  temperature 

363  °C 

j  AHC,  heat  of  combustion 

26000  kJ/kg 

AHV,  heat  of  vaporization 

1 108  kJ/kg 

Table  2  -  PMMA  Thermal  Properties  Used  in  the  FDS  Simulation 


5.2  A  Close  Look  at  PMMA’s  Thermal  Properties 

Poly(methyl  methacrylate)  (Figure  12),  which  is  usually  referred  to  as  PMMA, 
is  commonly  called  by  its  trade  name  PLEXIGLAS™.  As  a  polymer  underging  a 
polymerization  process,  PMMA  contains  no  cross-linking  between  long  chain 
molecules,  therefore  it  can  be  considered  as  a  thermoplastic.  PMMA  is  a  vinyl  polymer 
composed  of  vinyl  monomers  which  is  containing  carbon-carbon  double  bonds.  PMMA 
is  thus  the  product  of  the  polymerization  of  free  vinyl  radicals  called  monomer  methyl 
methacrylate. 
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Figure  12  -  PMMA 

In  a  recent  research  done  in  the  University  of  Maryland,  College  Park,  Steinhaus 

o  c 

and  Torero  evaluated  the  thermophysical  properties  of  PMMA.  The  research 
investigated  the  temperature  dependencies  of  many  thermal  properties  of  PMMA, 
which  is  determined  by  using  the  time  to  ignition  tig  and  the  time  to  pyrolysis  tp  as 
obtained  from  the  Forced  Flow  and  Flame  Spread  Test.  The  following  are  the  results 
presented  by  Steinhaus. 

Density  p 

It  is  known  that  the  density  decreases  linearly  with  increasing  temperature  and 
that  the  slope  is  constant.  The  following  expression  was  obtained  for  the  temperature 
dependent  density,  where  the  temperature  T  is  given  in  degree  Celsius: 

p  =  —7.3 1 6  •  1 0  4  •7’  + 1.4045  Equation  68 

cm  J 

As  denoted  by  the  above  expression,  there  is  approximate  20%  difference  for  the 
density  at  between  20  °C  and  380  °C. 

Specific  Heat,  cp 

The  specific  heat  has  to  be  broken  down  into  two  time  dependent  equations,  as 
dictated  by  the  glass  transition  phenomena  of  PMMA,  as  shown  in  Figure  1 3. 
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Figure  13-  Heat  Capacity  of  various  solid  and  molten  PMMA  35 


The  temperature  in  the  following  equations  are  given  in  Kelvin: 


c  =  2.04 1  •  1 0 2  -T~2  +  4.341  •  10“3  -r  +  0.173 


J 


gK 


Equation  69 


c  =  2.547  •10'-’  T  +  1.212 


J 


gK 


Equation  70 


Thermal  Conductivity,  k 

Many  investigators  have  conducted  research  on  the  thermal  conductivity  of 
PMMA,  and  their  results  are  summarized  in  Figure  14: 
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Figure  14-  Thermal  Conductivity  of  Various  Solid  and  Molten  PMMA  35 


The  obtained  curves  for  the  thermal  conductivity  as  a  function  of  time  for  the 
area  below  the  glass  transition  temperature  are  presented  in  Equation  71,  while  the  area 
above  the  glass  transition  temperature  is  expressed  in  Equation  72.  The  temperature  is 
in  °C. 


k  =  4.954-10'5  -T  +  0.1959 

k  = -2.318- 10‘3  -  T  +  0.2249 


W  ' 
m-°C 

’  W 
m-°C 


Equation  71 


Equation  72 


Pyrolysis  Temperature,  Tp 

A  series  of  tests  were  conducted  and  the  results  are  presented  in  Figure  15, 
where  the  literature  pyrolysis  temperature  is  given  by  the  manufacture  of  PMMA  (Cyro 
Industries)  used  by  Steinhaus35: 
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Figure  1 5  -  Thermocouple  Reading  on  the  Surface  of  the  Samples  at  the  onset  of 

pyrolysis  35 


Obviously  the  PMMA  used  by  Steinhaus  is  different  from  the  one  used  by  the 
FMRC  benchmark  experiment,  however,  Steinhaus'  study  showed  that  PMMA's 
thermal  properties  could  vary  by  a  great  percentage  between  ambient  temperature  and 
pyrolysis  temperature.  Such  nature  of  PMMA  causes  a  major  challenge  to  the  numerical 
simulation  of  fire  using  FDS,  which  uses  constant  thermal  properties  for  whole  course 
of  simulation,  regardless  of  the  sample  temperature 
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5.3  FMRC  Experimental  Results 


The  measured  total  heat  flux  histories  at  various  locations  on  the  PMMA  wall 
provide  an  important  insight  into  the  heat  flux  distributions  in  an  upward  propagating 
flame.  As  shown  in  Figure  16,  the  total  heat  flux  histories  along  various  heights  were 
given.  The  pyrolysis  height  reached  the  top  of  the  PMMA  wall  at  about  1200  seconds. 

The  heat  flux  distributions  are  clearly  going  through  three  phases.  Initially,  the 
flame  starts  with  a  triangle-like  profile.  As  the  flame  races  up  the  PMMA  wall,  the 
profiles  have  a  top-hat  distribution  with  peak  values  between  about  30-40  kW/m  . 
Finally,  after  the  pyrolysis  zone  reached  the  top  of  the  PMMA  wall,  the  profile  evolves 
toward  the  steady  state  shape.  Summing  up  the  total  energy  gives  a  radiation  fraction  of 
about  0.26. 

The  pyrolysis  height  history  is  given  in  Figure  17,  while  the  chemical  heat 
release  rate  history  is  plotted  in  Figure  18. 
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Figure  17  -  Pyrolysis  Height  History 
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Figure  18  -  Heat  Release  Rate  History 
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6  EVALUATION  OF  THE  FDS  PREDICTION  ON  THE  PMMA 


WALL  TEST 

6.1  FDS  Implementation  Issues 

The  FDS  application  of  the  existing  burning  rate  model  (Equation  34  ”  Equation 
37)  is  examined  hereinafter.  The  FDS  input  file  used  for  the  simulations  follows  in 
Figure  19.  The  FDS  simulation  domain  is  constructed  as  close  to  the  real  fire  test  as 
possible.  All  items  in  the  real  fire  test  is  constructed  into  the  Simulation  Domain, 
including  the  PMMA  vertical  fuel  wall,  the  marinite  panels  extended  from  the  fuel  wall 
and  the  steel  sheet  perpendicular  to  the  marinite  panels  which  is  used  to  minimize  the 
effects  of  room  drafts. 

According  to  Ma’s  finding  31  on  the  dependency  between  the  resolution  of  FDS 
domain  and  accuracy  of  FDS  prediction,  optimal  results  could  be  obtained  if  grid  cell 
was  set  at  about  5%  of  z*}  where  z*  is  shown  as  a  sole  function  of  the  fire  intensity: 

(  '  \a" 

2*  =  - — — j=  =  (q'd  J  5  ■  D  Equation  73 

In  the  FDS  simulation,  the  fire  intensity  ranges  from  about  20  kW  at  the 
beginning  as  a  result  of  the  ignition  device,  to  approximately  700  kW  at  the  stage  of 
quasi-steady  burning  over  the  whole  PMMA  wall.  The  corresponding  z*  ranges  from 
0.21  m  to  0.86  m,  while,  5%  of  z*  gives  a  value  between  1.0  cm  and  4.3  cm.  If  the 
resolution  of  the  FDS  domain  is  chosen  as  1.0  cm,  however,  there  will  be  3.6  millions 
cells  in  the  domain,  which  is  out  of  the  existing  FDS’s  handling  capacity.  Even  with  the 
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unofficial  working  version  of  FDS  V2.2  that  has  been  doubled  the  stacking  size,  such 
resolution  may  take  months  for  a  single  simulation  of  1500  seconds’  flame  spread  test. 

Then,  the  resolution  had  to  be  chosen  as  3.75  cm,  as  a  compromised  result 
between  fine  resolution  and  length  of  computer  time.  With  this  resolution,  the  total 
number  of  cells  is  approximately  66,000.  The  simulations  were  calculated  on  a  1.7  GHz 
Pentium  IV®  computer  with  2048  MB  of  RDRAM  memory.  It  took  this  computer 
approximately  150  hours  of  computer  time  for  a  simulation  of  1800  seconds’  flame 
spread. 

In  the  Experiment,  an  ignition  device  was  used  and  data  had  been  recorded  since 
the  ignition  device  was  turned  on.  The  ignition  device  was  a  long  hut  narrow  pan  of 
alcohol  pool  fire  along  the  wall  base,  ignited  by  a  match,  and  lasted  until  self-sustained 
flame  spread  could  be  achieved  on  the  wall.  If  to  perfectly  mimic  such  ignition  device 
with  az*  ”  10  mm,  an  unrealistic  ultrahigh  resolution  would  be  needed  and  the  number 
of  cells  would  be  as  high  as  10  millions.  Hence,  the  exact  simulation  of  the  ignition 
device  is  impossible. 

In  the  FDS  Simulation,  a  small  block  of  horizontal  propane  gas  burner  was 
placed  at  the  bottom  of  the  fuel  wall,  with  prescribed  initial  fire  size,  to  mimic  the  actual 
alcohol  ignitor.  The  fire  intensity  of  this  mimicked  ignitor  was  set  as  minimum  as 
needed  for  ignition  of  the  PMMA  wall,  and  this  was  achieved  by  hand  calculations  and 
trials.  The  resulted  fire  intensity  is  around  20  kW,  while  compared  to  the  actual  fire  test 
of  approximate  10  kW,  it  is  considered  acceptable  under  the  current  computational 
power  available  in  the  market. 
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FDS  INPUT  FILE 

&HEAD  CHID= ' omma '• TITLE- ' PMMA  Wall  Fire?1  / 

LGRID  IBAR-32,  JBAR-16,  KBAR-128  / 

tPDIM  XBAP-li|f  YBAR-0.60,  2BAR-5  / 

STIME  TWFIN-  3000./ 

LMISC  N FRAMES- 30 0 0 ,  DTCORE-30,  RESTART_FILE=  '  pmma  .  restart '  ,  REACTION-  1  MMA  '  / 
LREAC  ID- 1 MMA ' , 

FY I -  1 MMA  monomer,  C_5  H_8  0_2 1 

EPUM02=13542 

MW  FUEL-100 

NU  02-6 

NU^H20=4 

NU_  C02-5 

SOOT— YIELD— 0 . 022  / 

&SURF  ID-'PMMA’ 

HEAT_OF_VAPORIZATION=11D8 

HEAT__OF_COMBUSTION=260QO 

DELTA-0 . 025 

KS-0 .209 

ALPHA- 1 . 2E-7 

DENSITY-1190 

TMPIGN- 3  63  / 

&  SURF  I D- ! FIRE  1  ,  HRRPUA-2000  / 

SSURF  ID-  '  MARIN  I  TE  ' 

ALPHA- 1 . 14E-7 

KS-0. 12 

DELTA- 0.0254  / 

&5URF  1 D-1 SHEET  METAL' 

FYI-'IR  guaqe  Sheet  metal' 

C_DELTA_RHO- 4 . 7 

DELTA-0.0013  / 

SVENT  CB-  1  SEAR’  ,  S'JRF_1D-  1  OPEN  1  / 

SVENT  CB-  1  YBAR  1  ,  SURF__I  D=  ’  OPEN  1  / 

SVENT  xB- 0.0  0,0. 30, 0,0, 0.0, 0.00, 5, 00,  SURF_ID-  '  MAR1NITE  '  / 

SVENT  XB-0 . 90,  1  .  20,  0 . 0,  0  .  0,  0 . 00, 3 . 00,  SURF__ID=  '  MARIN  ITE  1  / 

6 VENT  cB= 'XBARO ' ,  SURF_ID= ' SHEET  METAL '  / 

SVENT  cB-'XBAR',  SURF_I D= ' SHEET  METAL'  / 

LVENT  XB=0 . 30, 0 . 90,  0 . 0,  0 . 03,  0 . 0,  0 . 0,  Sl)RF_I  D=  '  FIRE  '  / 

SVENT  XB-0. 30, 0 . 90, 0 . 0, 0 . 0, 0 . 0, 5 . 0,  SURF_ID= ' PMMA '  / 

SSLCF  XB=0 . 00, 1 . 20, 0 . 05, 0 . 05, 0 . 00, 5 . 00,  QUANTITY-; '  TEMPERATURE  '  / 

SSLCF  XB=0 . 60, 0 . 60, -0 . 05, 0 . 60, 0 . 00, 5 . 00,  QUANTITY- ' TEMPERATURE ’  / 

SBNDF  QUANTITY- 'WALL-TEMPERATURE'  / 

SBNDF  QUANTITY-  ’ BURNING^RATE ‘  / 

SBNDF  QUANTITY- '  KEA  T_FLUX’ 1  / 

SBNDF  QUANTITY-" '  CONVECTI VE_FLUX '  / 

Figure  19-  Input  File  for  FDS  Simulation  of  the  FM  Wall  Fire  Test 
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6.2  Comparison  of  the  Heat  Release  Rate  History 


Here,  the  experimental  Heat  Release  Rate  history  (Figure  18)  is  used  as  the 
ground  for  synchronizing.  By  matching  the  initial  HRR  of  34  kW  to  Figure  18,  the 
synchronizing  point  should  be  350th  seconds.  When  calculated  ignition  output  was  34 
kW,  this  was  “synchronized”  with  the  experiment  time  of  350  s.  After  synchronized, 
the  Heat  Release  Rate  history  comparison  is  shown  in  Figure  20.  Note  that  in  all  the 
following  comparison,  the  time  scales  are  synchronized  by  matching  the  350  seconds. 
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Figure  20  -  Comparison  of  HRR  History,  After  Synchronized  (34  kW) 


The  FDS  predicted  HRR  shown  in  Figure  20  are  obtained  by  the  equation 
Q  =  A//c  ‘  -  A cel/  •  The  FDS  prediction  is  in  excellent  agreement  with  the 

experimental  measurements 
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By  default,  the  FDS  uses  q’”=  AH 0  •  m”’)  (Equation  23)  for  the  Heat  Release 
Rate  per  unit  volume,  and  Q  =  'Vce!!  ■  However,  if  Equation  23  were  to  be  used, 

the  HRR  would  be  underestimated  by  S0%  maximum  (not  shown  in  Figure  20).  The 
reason  for  such  an  underestimated  heat  release  rate  is  that  there  is  a  limitation  with  the 
FDS  domain  described  in  Section  6.1.  The  actual  experimental  domain,  as  described  in 
Figure  10  and  Figure  11,  includes  a  3-m-high  extension  board  above  the  5-m-high 
PMMA  wall.  The  fire  product  collector  was  installed  above  the  extension  board,  and 
was  taking  account  the  whole  combustion  products.  However,  due  to  the  constrains  of 
time  frame  and  computational  power,  the  FDS  domain  used  in  this  work  is  S  m  high, 
just  including  the  PMMA  wall.  Thus,  the  combustion  occurred  above  the  S  m  PMMA 
wall  that  is  out  of  this  FDS  domain,  is  not  taken  into  account,  and  the  heat  released 
above  the  domain  is  not  taken  into  account  neither. 


6.3  Comparison  of  the  Flame  Spread 

Pyrolysis  height  predicted  by  the  FDS  is  observed  by  mass  flux  a0,  or  in  other 
word,  at  the  time  of  the  pyrolyzed  fuel  starting  to  inject  into  the  air.  The  result  for 
pyrolysis  height  is  shown  in  Figure  21.  The  code  successfully  simulated  the  upward 
flame  spread  in  a  satisfactory  way. 

The  surface  temperature  profile  predicted  by  FDS  is  presented  in  Figure  22.  It 
shows  the  heating  process  of  the  PMMA  wall  by  the  wall  fire.  And  it  conforms  to  the 
upward  flame  spread  phenomena,  because  the  surface  temperature  increases  when  time 
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and  height  increases.  The  maximum  value  of  the  surface  temperature  is  the  prescribed 
pyrolysis  temperature,  in  this  case,  363  °C. 
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Figure  21  -  Comparison  of  Flame  Spread 


400  i -  -  - 

1100s 


0.0  1.0  2.0  3.0  4.0  5.0 

Height  (m) 

Figure  22  -  PMMA  Surface  Temperature  Calculated  by  FDS 
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6.4  Comparison  of  the  Total  Heat  Flux 

At  various  height,  the  measured  maximum  values  for  total  heat  fluxes  are  within 
the  range  of  30-40  kW/m_.  In  FDS,  the  predictions  gave  the  maximum  values  as  high 
as  30  kW/m_.  However,  the  heat  flux  prediction  decreases  too  quickly  as  the  height 
increases.  The  finding  is  in  good  agreement  with  the  finding  in  the  no-spread  case 
described  in  Section  4.  The  comparison  of  total  heat  flux  distributions  at  various  times 
is  shown  in  Figure  23. 

The  FDS  predictions  on  heat  flux  histories  are  time-averaged  values  of  those 
values  of  the  nearest  ±  3  seconds,  in  order  to  minimize  the  moment  fluctuations. 
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Figure  23  -  Comparison  of  Total  Heat  Flux  Distribution  at  Various  Times 
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6.5  Comparison  of  the  Burning  Rate 

Predictions  of  burning  rate  is  taken  from  the  quasi-steady  stage  during  which  the 
whole  PMMA  wall  is  involved,  i.e.  1200s  “  1500s.  Although  Wu's  paper  33  didn't 
provide  experimental  data  on  the  burning  rate,  the  correlations  by  Ahmad  and  Faeth  15 
on  turbulent  convective  burning  of  vertical  plate  is  used  as  the  benchmark.  Orloff  and 
co-workers  at  FMRC  had  performed  a  smaller  flame-spread  experiment  on  the 
PMMA  wall  up  to  1  m  high  as  an  earlier  work  prior  to  Wu's.  Those  experimental 
measurements  are  also  served  as  benchmark  data  to  validate  the  FDS  predictions,  which 
is  time-averaged  for  the  data  between  1200s- 1500s.  The  comparison  is  presented  at 
Figure  24. 

Ahmad  and  Faeth's  correlation,  as  rewritten  by  Quintiere  37,  for  the  average 
burning  rate  m",  for  a  distance  x  measured  from  the  start  of  the  plate  is  given  by  the 
formula  below: 

—  ’*  'Cp  =  0.0285^/V*04  )>r073 /z  Equation  74 

Where: 


Equation  75 
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Equation  79 


Equation  80 


The  exact  thermal  properties  of  the  PMMA  as  shown  in  Table  2  are  used  in 
Ahmad  &  Faeth's  correlations,  plus  the  environment  variables  for  air  at  ambient 
temperature  as  shown  in  Table  3.  The  mass  burning  rate  (riiy )  from  Ahmad  &  Faeth's 
correlation,  as  shown  in  Figure  24,  are  discretized  values  from  the  averaged  mass 
burning  rate  ( m"F )  directly  given  by  Equation  80 


Thermal 

Properties 

Value  and  Unit 

K 

25E-3  W/m-s 

V 

15E-6m2/s 

Pr 

0.70 

P 

1.1  kg/m3 

cp 

1.05  kJ/kg-K 

T„ 

293  K 

Table  3  -  Air  Thermal  Properties  Used  in  Ahmad  &  Faeth's  Correlations 


Below,  Figure  24  is  the  results  when  the  whole  PMMA  wall  is  involved,  for  the 
quasi-steady  burning  stage.  The  results  are  time-averaged  for  the  period  between 
1500- 1800  s. 
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Figure  24  -  Comparison  of  Mass  Burning  Rate 


Again,  the  findings  from  Figure  24  shows  an  over-prediction  of  mass  burning 
rate  at  the  base  of  the  vertical  plate  and  somewhat  under-prediction  at  the  top  of  the 
vertical  plate.  This  finding  could  be  a  result  of  heat  flux  and  gas  temperature  prediction 
by  FDS  in  the  no-flame-spread  case  as  examined  in  Figure  6,  Figure  7  and  Figure  8, 
where  approximately  30%  over-prediction  of  gas  temperature  within  the  flame  region  is 
observed.  The  mass  burning  rate  for  the  region  between  4  m  and  5  m  is  underestimated. 
This  may  due  to  the  limitation  of  the  FDS  domain  used  in  this  simulation,  which  omits 
the  3-m  extension  board.  The  heat  flux  feedback  that  is  from  the  combustion  gas  out  of 
the  FDS  domain  is  not  taken  into  account,  making  a  substantially  underestimated  mass 
burning  rate 
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7  EVALUATION  OF  THE  FDS  PREDICTION  ON  THE  SMALL- 


SCALE  PMMA  WALL  TEST 

7.1  Small-scale  Flame  Spread  Experimental  Setup 

Two-dimensional  upward  flame  spread  and  subsequently  steady  turbulent 
burning  of  a  thermally  thick  vertical  fuel  surface  is  examined  experimentally  36.  Upward 
fire  spread  measurements  were  made  using  4.5  cm  thick,  41  cm  wide,  157  cm  high 
vertical  slabs  of  PMMA  cut  from  a  single  original  slab.  Since  the  upward  spread 
process  accelerates  and  is  therefore  intrinsically  unstable,  great  care  was  taken  to  insure 
truly  two-dimensional  spread  all  the  way  up  to  the  top  of  the  slab.  Sidewalls  (15  cm 
deep)  were  used  to  prevent  lateral  air  entrainment.  It  was  necessary  to  water-cool  the 
side-walls  to  prevent  excessive  radiant  heat  transfer  to  the  slab.  A  uniform  line  ignition 
was  achieved  at  the  bottom  edge  by  using  a  narrow  acetone  ignitor  tray.  The  averaged 
pyrolysis  height  (is.,  the  advancing  fuel  vaporization  front)  was  measured  by  visual 
observations  of  surface  bubbling  by  peering  through  the  rear  of  the  slab.  These  results 
were  reproduced  in  Figure  29. 


7.2  SOFIE  and  FDS  Calculations 

o 

Lewis,  Rubini  and  Moss  performed  a  field  modeling  of  flame  spread.  Using 
the  SOFIE  code  9,  a  simplified  model  of  non-charring  solid  pyrolysis  has  been  coupled 
with  detailed  descriptions  of  turbulent  combustion,  soot  production  and  radiative 
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exchange  to  create  a  field  model  for  fire  spread  prediction.  The  predictions  is 


demonstrated  by  comparison  with  established  two-dimensional  experiments  on  PMMA 
slabs  reported  in  the  literature  by  Orloff  et  al  36.  The  SOFIE  Code  uses  a  two-equation 
k-E  turbulence  model  for  modeling  of  buoyant  turbulent  flow,  while  the  FDS  uses  the 
LES  as  its  turbulence  model.  It  was  reported  that  the  level  of  agreement  between 
measurements  of  flame  spread  velocity,  mass  loss  rate  and  incident  heat  flux  is 
generally  good. 

In  order  to  perform  a  cross-comparison  for  the  predictions  by  FDS  and  SOFIE 
and  the  experimental  measurements,  a  FDS  domain  is  designed  and  the  simulation 
method  is  similar  to  the  large-scale,  5-m-high  flame  spread  case.  The  only  exception  is 
that  the  thermal  properties  that  is  used  in  this  small-scale  domain  is  chosen  to  be  as 
close  as  possible  with  those  used  by  SOFIE  calculations  8. 

A  baseline  computation  has  been  defined  for  model  comparison  purposes  on  a 
117  x  32  mesh  (H  x  W)  on  the  PMMA  surface.  Atop  of  the  PMMA  wall  is  a  1-m-higb 
gypsum  board  extension,  served  as  the  flame  holder.  The  FDS  domain  consists  of  32  x 
x  12  x  256  by  its  width,  depth  and  height,  and  there  were  over  75,000  cells  in  the 
calculation  domain.  The  resolution  of  the  domain  is  as  small  as  1,3  cm,  and  falls  into 
the  range  of  5%  of  z*,  the  characteristic  length  of  the  firepower. 


7.3  The  PMMA  Thermal  Properties 

In  FDS,  fixed  values  are  used  for  thermal  properties  to  be  used  as  inputs  for  the 
calculation.  However,  it  seems  that  SOFIE  uses  the  variable  thermal  properties  for  the 
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PMMA  thermal  properties  38.  Then,  averaged  values  for  k  and  a  are  picked  for  the  use 
of  FDS  calculations.  The  thermal  properties  are  shown  in  Table  4,  and  the  comparison 
of  values  of  thermal  conductivity  and  thermal  diffusivity  are  shown  in  Figure  25  and 
Figure  26  respectively. 


Thermal  properties 

SOFIE 38 

FDS 

Specific  heat 
(Cp) 

2.374  T+  1100 
J/Kg/K 

- 

Density 

(p> 

1190  Kg/m3 

1190  Kg/m3 

Conductivity 

m 

2.49e-4  T  +  0.1 18 
W/m/K 

0.23366  W/m-K 

Thermal  diffusivity 
(a) 

- 

8.9141  E-8  m2/s 

Heat  of  vaporization 

(  AHV) 

1.108  MJ/kg 

l,108MJ/kg 

Ignition  temperature 
(TP) 

636  K 

636  K 

Table  4  -  Thermal  Properties  Used  In  SOFIE  And  FDS 
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Thermal  Diffusivity  (m  2/s)  Thermal  Conductivity,  k  (J/kg-K) 
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Figure  25  -  Comparison  Of  k 
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Figure  26  -  Comparison  Of  a 
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7.4  Comparison  of  the  Results 


The  FDS  code  (V  2.0)  does  not  finish  the  calculations,  as  it  ends  automatically, 
but  accidentally,  when  the  HRR  increases  in  a  very  rapid  rate.  In  all  the  cases,  the 
calculations  ended  with  unknown  reasons,  regardless  of  the  thermal  properties  and 
domain  resolutions.  In  the  FDS  calculations,  various  thermal  properties  are  used  one  by 
one  for  numerous  trials:  reduced  k  (0.209  W/m-K);  increased  a  (1.2E-7  m2/s);  increased 
AHy  (1.6  MJ/kg),  etc.  In  addition,  the  resolution  has  been  modified  and  coarse 
resolution  (cell  number  reduced  50%)  has  been  tried.  In  all  these  calculations,  the 
similar  “blow  up”  phenomena  are  occurred  and  FDS  ends  the  calculations  accidentally. 
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Figure  27  -  Heat  Release  Rate  Given  By  FDS  Calculation  For  The  Small-Scale  Test 
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Figure  28  -  Cross  Comparison  For  Heat  Release  Rate  Given  By  FDS  Calculation  For 

The  Large-  And  Small-Scale  Tests 


In  Figure  28,  the  FDS-calculated  results  for  large  and  small-scale  flame  spread 
tests  are  compared  with  the  5-m  experiment.  There  are  differences  for  the  geometry 
settings  between  the  large  and  small-scale  tests.  One  of  those  is  that  between  the 
perpendicular  side-plates,  the  large-scale  test  installed  two  60-cm-wide  gypsum  boards 
beside  the  PMMA  wall,  while  the  small-scale  test  didn't.  This  may  cause  that  a  higher 
heat  flux  on  the  fuel  surface  for  the  small-scale  test  and  therefore  faster  flame-spread. 

The  FDS  domain  for  the  large-scale  test,  due  to  the  constrain  of  long  computing 
time,  uses  a  coarse  resolution  (3.7  cm),  and  causes  inconsistencies  in  simulating  the 
ignition  device  which  is  located  at  the  base  of  the  PMMA  wall.  To  address  this 
problem,  the  FDS-calculated  results  are  synchronized  to  match  the  same  heat  release 
rate  of  the  actual  ignition  device.  When  calculated  ignition  output  was  34  kW,  this  was 
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“synchronized  with  the  experiment  time  of  350  s.  However,  the  FDS-domain  for  the 
small-scale  test  uses  a  much  better  resolution  (1.3  cm),  so  that  the  ignition  device  can  be 
simulated  rather  accurately.  In  Figure  28,  the  FDS  result  for  large-scale  test  is  the 
synchronized  one,  while  the  FDS  result  for  the  small-scale  test  is  its  original  result. 
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Figure  29  -  Cross  Comparison  Between  Predicted  Flame  Spread  Velocity  By  SOFIE 
And  FDS  And  Experimental  Measurements  By  Orloff 


Figure  29  shows  the  FDS-predicted  flame  spread  velocity  is  higher  than  the 
experimental  results.  The  flame  spread  velocity,  v,  is  related  to  the  surface  temperature 

(val /(Tjg  -TJ2).  When  T  — *  Tlg ,  v  -*■  00  .  Accelerated  flame  spread  between  490s 

and  542s  is  shown  in  Figure  29.  The  relationship  of  Ts  and  Tjg  is  also  shown  in  Figure 
30. 
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Figure  30  -  FDS-Calculated  Surface  Temperature  Vs.  Height,  Varied  By  Time 

Compared  to  Figure  22,  the  surface  temperature  calculated  by  FDS  appears  to  be 
much  faster  increase.  This  may  be  the  result  of  higher  heat  flux  as  shown  in  Figure  31 
and  Figure  32. 

In  Figure  3  1  and  Figure  32,  the  steady  stage  for  Orloffs  experimental  results 
and  SOFIE-calculated  results  are  time-averaged  values  at  an  “after  spread”  stage,  when 
the  entire  PMMA  wall  is  burned  and  is  maintained  quasi-steadily.  These  results  are 
compared  to  the  FDS-calculated  results  that  are  taken  for  developing  times,  and  the 
FDS  simulation  never  reaches  a  steady  stage  because  it  stops  unexpectedly. 
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Figure  3  1  —  Incident  Convective  Heat  Flux  Comparison 

Figure  31  shows  the  results  for  convective  heat  flux,  which  is  not  dominant 
comparing  to  radiative  heat  flux.  The  Orloffs  experimental  results  give  that  the 
convective  heat  flux,  in  general,  consists  of  less  than  20%  of  the  total  incident  heat  flux. 
The  FDS  gives  satisfactory  prediction  for  the  convective  heat  flux. 
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Figure  32  -  Incident  Radiative  Heat  Flux  Comparison  Between  Predictions  And  The 

Experimental  Values 

Figure  32  shows  the  results  for  radiative  heat  flux.  The  trend  of  FDS  is  very 
different  from  the  one  given  by  SOFIE  or  the  experimental  result.  The  radiative  heat 
flux,  q"r,  is  related  to  the  flame  temperature  7/  in:  q”r  oc  EoTf4 .  As  the  flame 
temperature  is  expected  to  be  approximately  constant,  and  the  emissivity 
e  «  1-  w  where  6  is  the  boundary  layer  thickness  which  increases  with  height. 
Hence,  the  radiative  heat  flux  should  increase  as  the  height  increases.  Figure  32  shows 
SOFIE  gives  the  proper  results  while  FDS  does  not.  The  reason  might  be  that  FDS  is 
giving  an  incorrect  flame  temperature,  as  shown  in  Figure  6 
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8  SUMMARY  AND  CONCLUSIONS 


Flame  spread  and  burning  rate  is  examined  in  a  5-m  high  PMMA  vertical  slab, 
using  a  CFD  code,  the  Fire  Dynamics  Simulator  (FDS).  Coupled  descriptions  of  the 
controlling  mechanisms  of  hydrodynamics  model  for  buoyancy-driven  flow,  mixture- 
fraction  model  for  combustion,  convection  model,  banded  gray-gases  transportation 
model  for  radiation,  one-dimensional  conduction  model  for  pyrolysis  and  burning  rate 
model  are  included.  The  approaches  for  pyrolysis  and  burning  rate  model  are  evaluated 
in  sub-model  level,  and  the  refined  algorithms  are  presented.  The  new  algorithms  bring 
in  the  concept  of  thermal  penetration  depth.  One  of  the  major  benefits  of  this  new  term 
is  that  the  thermal  penetration  depth  is  independent  of  the  resolution  of  FDS  domain, 
thus  reducing  the  dependency  of  FDS  prediction  on  surface  temperature  and  burning 
rate.  The  proposed  algorithm  is  not  incorporated  into  the  FDS  code  and  is  a  subject  of 
future  study  for  quantifying  its  computational  benefit,  and  the  validation  works  are 
based  on  the  existing  models. 

The  models  are  first  tested  on  a  non-spread  case  to  examine  the  accuracy  of 
prediction  for  heat  flux  and  gas  temperature  distribution.  Results  were  compared  with 
the  literature  and  the  radiation  sub-model  is  tested.  The  FDS/LES  model  gives  good 
results  for  flame  height.  The  results  for  temperature,  wall  heat  flux,  and  flame  height  are 
in  qualitative  agreement  with  data.  The  temperature  and  the  heat  flux  are  within  about 
30%  of  measured  results  for  peak  values  near  the  base  of  the  wall,  but  both  decay  much 
quicker  than  the  experimental  results. 
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Further  predictions  are  made  on  a  more  complex  large-scale  flame  spread  case 
involved  a  PMMA  wall,  which  is  5  m  tall  and  0.6  m  wide,  with  perpendicular  steel  draft 
containing  the  flame  sheets.  Predictions  on  heat  flux,  gas  and  solid  temperature,  flame 
spread,  and  mass  burning  rate  are  evaluated.  The  FDS  model  gives  satisfactory  flame 
spread  predictions.  The  heat  release  rates  are  within  50%  of  the  measured  experimental 
values.  The  total  heat  flux  predictions  are  within  30%  of  measured  results  at  the  lower 
part  of  the  PMMA  wall  but  it  drops  quickly,  and  the  deviation  increases  as  the  height 
increases.  This  is  in  agreement  with  the  indications  of  gas  temperature  and  heat  flux  in 
the  non-spread  case.  The  model  over-predicts  the  burning  rate  at  the  lower  part  of  the 
PMMA  wall.  But  the  model  gives  an  excellent  result  for  the  2-4  m  region,  as  the 
predictions  are  within  10%  range  of  Ahmad  and  Faeth’s  correlations.  The  poor 
predictions  for  the  top  region  (4-5  m)  are  probably  due  to  the  reason  that  the 
computational  domain  lacks  the  upper  extended  steel  sheets  (as  shown  in  Figure  10) 
that  act  as  a  “holder”  to  reduce  flame  fluctuations,  while  the  experimental  setup  did 
equip  such  extensions. 

A  cross-comparison  is  made  for  FDS  and  SOFIE  for  a  small-scale  flame  spread 
test.  It  appears  that  SOFIE  does  a  better  job,  and  FDS  gives  inconsistent  results  for  the 
prediction  of  radiative  heat  flux.  It  seems  that  FDS  suffers  from  poor  accuracy  for  near 
fuel  base  where  fuel  is  rich.  Here  the  gas  temperature  is  overestimated,  and  the  data  is 
underestimated  when  the  system  becomes  fuel  lean  at  the  far  field. 

Generally,  reasonable  agreement  is  achieved  with  experiments  at  levels  of 
spatial  resolution  that  are  compatible  with  the  computational  demands  of  field  model 
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predictions  in  compartment  fire  scenarios.  The  strength  of  FDS  appears  in  that 
entrainment  is  calculated  accurately  since  flame  extent  is  well  predicted.  The  deficiency 
of  FDS  seems  to  be  that  combustion  details  are  distorted,  giving  higher  gas  temperature 
in  region  of  the  fuel-rich  area  and  lower  gas  temperature  for  the  fuel-lean  area. 
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APPENDIX  A 


Proposed  Algorithm  of  Predicting  Solid  Fuel's  Surface  Temperature 

The  numerical  solution  was  developed  and  provided  informally  by  Professor 
James  G.  Quintiere  from  the  Department  of  Fire  Protection  Engineering  at  the 
University  of  Maryland.  The  solution  was  also  given  in  the  dissertation  of  Francisco 
Joglar-Billoch  at  the  University  of  Maryland39.  Here  is  a  reproduction  of  the  solution. 

Recalling  the  solution  for  the  differential  equation  of  1-D  conduction  problem, 
the  following  solution  was  developed  to  estimate  the  surface  temperature  when  a  solid 
is  exposing  to  a  time-related  net  surface  heat  flux  q"net{ x)  =  q"c  +  q"r  ~q"RR  : 

1  1  ' "  (  } 

Ts  -  Ta  =  —7 - C^n.——dx  Equation  81 

■q  Jtkpc  \  t  —x 

Rewrite  as, 

t  •  if  (  \ 

AT  =  V  C^Mix  Equation  82 

JoV/-r 


where  the  variable  W  is  used  to  represent  the  constant  term  of  the  equation  1 /^jjikpc  . 
Equation  81  can  be  represented  into  a  "time  step"  form: 


y  /  (rt  + 1 ) 


Equation  83 


Using  the  trapezoidal  rule,  the  first  integral  in  the  right  hand  side  of  the  equation 
can  be  solved.  However,  the  second  integral  has  a  singularity  since  x  =  t  +  1 .  This 
singularity  can  be  removed  by  first  considering  Figure  33: 
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Figure  33  -  Numerical  Solution  of  Equation  83 


2udu  =  -dt 

At  t  =t(n),  u2  =t(n+l)  +  t(n)  or  u  =  VAt 

At  x  =  t(nTl),  u2  =  0 

Therefore,  in  the  transformed  variable  u, 
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0  'tt  /lrt  +  L  2  \  .  x 

-u  )du  or  /2=2L:J--m2>u 
£  u  0 


Equation  87 


and  using  again  the  trapezoidal  rule, 


|VaT  or 

Finally,  the  solution  for  Equation  81  is  given  by: 


Equation  88 


C  =7;  + 


=(/,  +  /2) 


hi  •  A'pc 


Equation  89 
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